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Abstract. We present and discuss new determinations of metallicity, rotation, age, kinematics, and Galactic orbits 
for a complete, magnitude-limited, and kinematically unbiased sample of 16,682 nearby F and G dwarf stars. 
Our ~63,000 new, accurate radial-velocity observations for nearly 13,500 stars allow identification of most of 
the binary stars in the sample and, together with published uvbyf3 photometry, Hipparcos parallaxes, Tycho-2 
proper motions, and a few earlier radial velocities, complete the kinematic information for 14,139 stars. These 
high-quality velocity data are supplemented by effective temperatures and metallicities newly derived from recent 
and/or revised calibrations. The remaining stars either lack Hipparcos data or have fast rotation. 
A major effort has been devoted to the determination of new isochrone ages for all stars for which this is possible. 
Particular attention has been given to a realistic treatment of statistical biases and error estimates, as standard 
techniques tend to underestimate these effects and introduce spurious features in the age distributions. Our ages 
agree well with those by Edvardsson et al. IJ1993J, despite several astrophysical and computational improvements 
since then. We demonstrate, however, how strong observational and theoretical biases cause the distribution of 
the observed ages to be very different from that of the true age distribution of the sample. 

Among the many basic relations of the Galactic disk that can be reinvestigated from the data presented here, we 
revisit the metallicity distribution of the G dwarfs and the age-metallicity, age-velocity, and metallicity-velocity 
relations of the Solar neighbourhood. Our first results confirm the lack of metal-poor G dwarfs relative to closed- 
box model predictions (the "G dwarf problem"), the existence of radial metallicity gradients in the disk, the small 
change in mean metallicity of the thin disk since its formation and the substantial scatter in metallicity at all 
ages, and the continuing kinematic heating of the thin disk with an efficiency consistent with that expected for a 
combination of spiral arms and giant molecular clouds. Distinct features in the distribution of the V component 
of the space motion are extended in age and metallicity, corresponding to the effects of stochastic spiral waves 
rather than classical moving groups, and may complicate the identification of thick-disk stars from kinematic 
criteria. More advanced analyses of this rich material will require careful simulations of the selection criteria for 
the sample and the distribution of observational errors. 

Key words. Galaxy: disk - Galaxy: solar neighbourhood - Galaxy: stellar content - Galaxy: kinematics and 
dynamics - Galaxy: evolution - Stars: fundamental parameters 



1. Background and motivation 

The Solar neighbourhood is the benchmark test for models 
of the Galactic disk: The stars in a sample volume around 
the Sun provide a first estimate of the mass density of the 
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Galactic disk near the plane. Their distribution in age is 
our record of the star formation history of the disk. Their 
overall and detailed heavy-element abundances as func- 
tions of age are the fossil record of the chemical evolution 
and enrichment history of the disk. Finally, the space mo- 
tions and Galactic orbits of the stars as functions of age 
are our clues to the parallel dynamical evolution of the 
Galaxy and the degree of mixing of stellar populations 
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from different regions of the disk (see, e.g. the recent re- 
view by Freeman & Bland-Hawthorn I2002f) . 

Providing the basic data for the stars of the 
Solar neighbourhood - ages, metallicities, velocities, and 
Galactic orbits - may seem the simplest observational 
task of all. Yet, the identification of the nearest stars 
and the data needed for a proper characterisation of their 
main astrophysical parameters remain seriously incom- 
plete. Moreover, stars have often been selected for obser- 
vation from criteria which induce correlations between the 
observed parameters that reflect the characteristics of the 
selection process more strongly than those of the Galaxy. 
Intrinsic parameters such as the frequency, age, and metal- 
licity of groups of stars are intertwined with such easily 
observable properties as brightness, colours, and proper 
motions in ways that can make it difficult or impossible 
to retrieve the former from the latter. Kinematic selection 
biases are particularly pernicious in this regard. 

F- and G-type dwarf stars are favourite tracer popula- 
tions of the history of the disk. They are relatively numer- 
ous and sufficiently long-lived to survive from the forma- 
tion of the disk; their convective atmospheres reflect their 
initial chemical composition; and ages can be estimated 
for at least the more evolved stars by comparison with 
stellar evolution models. Photometry in the Stromgren 
uvbyf3 system is an efficient means to derive their intrin- 
sic properties from observation (Stromgren 119631 IT587). 
Recent studies based mainly on uvby/3 photometry are, 
e.g. Feltzing et al. pi) . Holmberg l|2TM|l . and J0r gensen 
(2000), while more detailed spectroscopic studies of the 
chemical history of selected subsets of stars have been 
made by, e.g. Edvardsson et al. I|1993[> . Fuhrmann l|1998|l . 
Reddy et al. PTMjl . and Bensby et al. l|2lIU3|) . 

Extensive uvbyfi photometric surveys of the nearby F 
and G stars have been performed by Olsen <|1983l fl993 
I1994al I1994bj) . Accurate parallaxes and proper motions 
have become available for large numbers of these stars 
from the Hipparcos (ESA I1997[) and Tycho-2 (H0g et al. 
2000) catalogues. The bottleneck so far has been the cor- 
responding radial- velocity data. These are needed, first, to 
complete the three-dimensional space motions for all the 
stars and improve the statistical accuracy of the derived 
age-velocity relations. They also serve to identify those 
stars which just happen to pass near the Sun at present, 
but were formed elsewhere and have witnessed another 
evolutionary history than that of the Solar circle. Finally, 
and perhaps most importantly, repeated radial-velocity 
measurements allow identification of the large fraction of 
stars in the photometric samples which are binaries, and 
for which the derived astrophysical and kinematical data 
will be inaccurate and potentially misleading. 

The key contribution of this paper consists of new, ac- 
curate, radial velocities for an all-sky, magnitude-limited 
and kinematically unbiased sample of ^13,500 nearby F 
and G stars, based on ^63,000 individual photoelectric 
observations. Because this data set is unlikely to be super- 
seded in essential respects until the results of the GAIA 
mission (Perryman et al. 2001 ; and/or the RAVE project 



(Stcinmctz 2003) become available, we have also recal- 
ibrated and redetermined the astrophysical parameters 
(T e ff, M v , and [Fe/H]) for all stars in our sample. Much 
effort has been devoted to the fundamental issue of deter- 
mining reliable isochrone ages for as many stars as possi- 
ble, and we believe that a rather more realistic assessment 
of the errors of such ages has been obtained. Finally, we 
have computed individual Galactic orbits for all stars with 
adequate data. 

Table 1. Table 1 is only available in electronic form at 
the CDS. The first two pages of the table, listing the first 
100 stars, are given at the end of this paper as sample of 
its content and format. 



The resulting data set should place a wide range of 
studies of the evolution of the Galactic disk on a new 
and considerably improved footing. We caution, however, 
that no such thing as a fully unbiased sample exists in 
Galactic astronomy, and the reader is strongly urged to 
carefully study Sect. [SJ where we discuss the key issue of 
completeness of the data, especially as regards the derived 
ages and masses. The catalogue with the complete data set 
for 16,682 stars will be available electronically 1 (sample 
pages are shown in Table |2J). 

The complementarity between the present study and 
the much-quoted paper by Edvardsson et al. (1993 ) de- 
serves clarification. The 189 stars studied by Edvardsson 
et al. lfH?9Hl were drawn from the kinematically unbiased 
sample presented here, selecting equal numbers of stars 
in 10 metallicity bins through the range of [Fe/H] seen 
in the (thin and thick) disk. Evolved F dwarfs were se- 
lected, excluding known binaries and fast-rotating stars, 
so that interstellar reddening and isochrone ages could be 
determined. Thus, unevolved stars were avoided, young 
metal-poor and old metal-rich stars were excluded a pri- 
ori by the colour cutoffs used, and the sample was strongly 
biased in favour of metal-poor stars. The sample discussed 
here was explicitly designed to alleviate these selection bi- 
ases and further benefits from the advent of the Hipparcos 
and Tycho-2 data but, for obvious reasons, lacks the de- 
tailed, precise, and homogeneous spectroscopy that was 
the centrepiece of the study by Edvardsson et al. (1993 ). 

The present paper is organised as follows: Sect.Elclar- 
ifies the selection criteria used to define the sample, and 
Sect. |3 describes the new radial velocities and other basic 
observational material available for the stars. Sect. 0] ex- 
plains the (partly new) calibrations and analysis methods 
we have used to derive reliable astrophysical parameters 
from the raw data, and Sect. El discusses the complete- 
ness and various biases in the resulting parameter sets. 
We briefly rediscuss some of the classical diagnostic di- 

Tables 1 and 2 are available in electronic form by anony- 
mous ftp to the CDS at cdsarc.u-strasb.fr (130.79.128.5) or via 
|http://cdsweb.u-strasb.fr/cgi-bin/qcat?J/A+A/418/989| 
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agrams in Sect. [S]and finally summarise our conclusions 
and some of the prospects for the future in Sect.0 

2. Sample definition 

From the outset, it was clear that the definitive selection 
and characterisation of the nearby F and G dwarfs should 
be made from a comprehensive photometric survey in the 
Stromgren uvby/3 system. But in order to undertake such 
a survey, an observing list is needed. Thus, the first step in 
the project was to establish an apparent-magnitude lim- 
ited, kinematically unbiased all-sky sample of stars within 
which the F and G dwarfs would also be volume-complete 
to a sufficiently large distance (~40 pc). 

For this, we chose the HD catalogue (Cannon & 
Pickering 1918-21J, the only all-sky spectral catalogue 
then available. All A5-G0 stars brighter than m V i S = 8.3, 
GO stars in the interval 8.30< m V i S <8.40, and all G5 or 
just G stars (no subtypes) brighter than m vis = 8.6 were 
selected. On the one hand, the A5 spectral-type limit is 
early enough to include even quite metal-poor F stars; on 
the other hand, most KO-type stars in the HD catalogue 
are giant stars for which luminosities, distances, tangential 
velocities, and ages cannot be reliably determined. 

It was desirable to also include any old, metal-rich 
dwarfs and thus improve the observational basis for a re- 
assessment of the "G-dwarf problem" . However, observing 
the vast numbers of HD K-type giants is an inefficient way 
to identify such stars. We therefore added all 1,277 GOV- 
K2V stars south of S = -26° from the Michigan Spectral 
Catalogues (Houk & Cowley 1975, Houk 1978, 1982) that 
were previously unobserved, but estimated from the spec- 
tral types and m pg to be within a distance of 50 pc, 
with generous allowance for the uncertainties. This sam- 
ple should be complete to about the same distance limit 
(40 pc) as the rest of the G-dwarf sample, but covers only 
part (28%) of the entire sky. 

The vast majority of all these stars (^30,000 in to- 
tal) were then observed in the uvby/3 photometric surveys 
described below, and the sample for the radial-velocity 
programme was defined from the measured photometric 
indices. Subsequently, the Hipparcos and Tycho-2 cata- 
logues have provided accurate parallaxes and proper mo- 
tions for nearly all the stars for which radial velocities had 
been measured. 

2.1. Initial photometric survey 

From the observing lists assembled as described above, 
nearly all stars were observed at least once in the 
Stromg ren uvby/3 system by Olsen l(T3S3l ITMHl 11994a! 
I1994b|) . The resulting catalogues were merged with pre- 
vious sources of uvby/3 photometry (Stromgren & Perry 
IT5?)5l Crawford et al.lll^[ll^lll9Vlal[Wlblll9V2l[llI73l 
Gr0nbech & Olsen IT5?6l 115771 and Olsen & Perrv lT38ll) . 
The combined FG photometric catalogue thus constructed 
contains a total of 30,465 stars. 
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This all-sky sample is complete to the magnitude lim- 
its described above and should be volume complete for 
the F and G dwarfs to a distance of ~40 pc. In the south- 
ern cap where modern MK spectral types are available, 
the G0V-K2V stars within the same distance are also in- 
cluded. This homogeneous, complete, and kinematically 
unbiased database was used to select stars for the radial- 
velocity programme, using the criteria described below. 




Fig. 1. Distribution on the sky of the 16,682 programme 
stars. Top: The 14,139 stars with radial velocity data (note 
the overdensity of stars in the Hyades cluster and south 
of (5 = —26°). Bottom: The 2,543 stars with no radial 
velocity. 



2.2. Sample definition for the catalogue 

The sample of stars for which radial velocities have been 
obtained was defined from the complete FG catalogue. 
Slightly generous limits in photometry space have been 
adopted to allow a selection of FG dwarf stars based on 
physical criteria like mass, temperature, metallicity, etc. 

The following four samples were selected, merged, and 
cleaned of duplicate entries: 



4 



B. Nordstrom et al.: The Geneva-Copenhagen survey of the Solar neighbourhood 



1. All stars for which the F-star calibrations of Crawford 
(1975) and Olsen (1988) are valid. 

2. All stars with no (3 value and 0.240 <b-y< 0.460, 
[mi] > 0.120, 6ci < 0.400, and V < 9.600 

3. All stars with 0.205 < b -y < 0.240, [mi] > 0.120, 
6a < 0.400, and V < 9.600. 

4. All stars for which the G- and K-star calibrations of 
Olsen (1984) are valid. 

Criterion 3 ensures that no metal poor F-stars will 
be lost on the hot side of the F-type stars. Both crite- 
rion 2 and 3 compensate for some missing (3 observations. 
Criteria 1 and 4 and the V- limit 9.6 ensure that a number 
of fainter stars observed for calibration purposes are also 
included. After removal of a few known supergiants and 
other irrelevant objects, the list contains a total of 16,682 
objects. 

The distribution of the full sample of 16,682 stars over 
the sky is shown in Fig. ^ in equatorial coordinates and 
in an equi-area projection. Note the concentration of stars 
in the Hyades - observed with special care for calibration 
purposes - and the addition of the latest-type dwarfs south 
of 8 = —26°. Apart from these features, the sample is very 
uniformly distributed on the sky. 

3. Observational data 

For reference, we summarise in the following the basic 
observational data underlying the astrophysical and kine- 
matical parameters derived for the programme stars. Our 
new radial-velocity observations are described in detail; 
for other data, the relevant sources are given. 

3.1. Stromgren uvbyf3 photometry 

The merged catalogue of uvby[3 photometry from which 
our programme stars were selected was briefly described 
above (Sect. I2.2|l : much additional explanation and ex- 
tensive notes are given in the original catalogues (Olsen 
IT5£3l lEMl I1994al I1994b|l . With time, the (3 observations 
were extended to include all stars with a b-y value indicat- 
ing that (3 might be a useful reddening-free temperature 
indicator, except that (3 observations for a small fraction 
of the northern sky are still pending. The distribution of 
the sample in b-y colour is shown in Fig. 

3.2. Radial velocities 

When the first phase of this project was initiated, not 
even the Bright Star Catalog (Hoffleit & Jaschek IT§g2"jl 
had complete radial-velocity coverage. Of the ^1,500 
missing southern stars in that catalogue, the early-type 
half were observed with conventional spectrographic tech- 
niques (Andersen & Nordstrom I1983al I1983bl Nordstrom 
& Andersen HT)g5|) . the late- type half with CORAVEL 
(Andersen et al. 119851 see also Sect. I3.2TT|) . 

For the fainter stars, new observations were obtained 
as detailed below - altogether, a total of 62,993 new radial- 
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Fig. 2. Distribution of the whole sample in b-y colour. 

velocity observations of 13,464 programme stars. Adding 
earlier literature data, complete kinematical information 
is available for a total of 14,139 stars. 

3.2.1. CORAVEL observations 

The bulk of the radial-velocity data presented here was 
obtained with the photoelectric cross-correlation spec- 
trometers CORAVEL (Baranne et al. 113791 Mayor [T^SJ). 
Operated at the Swiss 1-m telescope at Observatoire de 
Haute-Provence, France, and the Danish 1.5-m telescope 
at ESO, La Silla, the two CORAVELs cover the entire sky 
between them, and their fixed, late-type cross-correlation 
template spectra efficiently match the spectra of the large 
majority of our programme stars. 

Initially, specific observing programmes were targeted 
to primarily cover the thick-disk stars which were assumed 
to be very old, and the evolved thin-disk stars for which 
ages can be determined. When a separate programme to 
observe all the late-type southern stars of the Hipparcos 
survey was initiated (Udry et al. [ESTf , most of the re- 
maining southern unevolved F stars were included as well. 
Subsequently, a good fraction of the northern half of the 
stars has been also observed in a separate Geneva pro- 
gramme on the Hipparcos stars. 

In all programmes, two or more observations were 
made for almost all stars over a substantial time base. 
This allows to define more reliable mean velocities, but 
also to identify most of the spectroscopic binaries which, 
if unrecognised, yield misleading astrophysical parameters 
from the observed magnitudes and colour indices. Between 
the two telescopes, a total of 60,476 CORAVEL observa- 
tions have been made of 12,941 of the programme stars 
discussed in this paper - some 1,000 nights' worth of data. 
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The catalogue presented here contains the mean radial 
velocity for each star together with the summary data 
on the observations as described in the Appendix. Our 
computations of the observational errors, criterion for de- 
tecting variable (i.e. binary) stars, and our treatment of 
double- lined binaries are discussed in Sect. 13. 2l4l below. 

Many stars on the main programme are primaries of 
close double stars. CORAVEL observations were made of 
the fainter companions to many of these stars in order 
to ascertain whether they are physically bound or merely 
optical pairs. These data will be made available separately 
and are not discussed further here. 

3.2.2. CfA observations 

The fixed-resolution CORAVEL mask is optimised for 
sharp-lined spectra, and the cross-correlation profile 
rapidly becomes too broad and shallow to yield accu- 
rate radial velocities for stars rotating faster than 40-50 
km s _1 , as is the case for a large fraction of stars with 
0.20 <b-y< 0.27. 

In order to recover as many as possible of the early F 
stars of the sample, several hundred stars rotating too 
rapidly for CORAVEL were observed with the digital 
spectrometers (Latham 1985) of the Harvard-Smithsonian 
Center for Astrophysics (CfA). These instruments yield 
accurate results for single stars with rotational velocities 
up to ~ 120 km s -1 (Nordstrom et al. 1994) and also per- 
form well on double-lined spectra with two-dimensional 
cross-correlation techniques (Latham et al. 1996). 

CfA radial velocities for 595 stars were published by 
Nordstrom et al. (1997b) and have been used in the 
present catalogue when CORAVEL data were missing or 
less accurate (2,517 observations of 523 stars). Rapidly ro- 
tating stars south of declination -40° cannot be reached by 
the CfA instruments and thus have no new radial- velocity 
data. 
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3.2.3. Literature data 

When this programme was initiated, radial-velocity data 
existed in the literature essentially only for stars in the 
Bright Star Catalog. As the earlier data were mostly of 
reasonably good quality, these bright stars have generally 
not been reobserved. Literature data for a total of 675 
such stars and others not covered by the new programmes 
have been taken, as far as possible, from the compilation 
of Barbier-Brossat & Figon (2000), and bring the total 
number of stars with radial velocity data to 14,139. 

3.2.4. Variability criterion and binary detection 

Each radial-velocity observation is associated with an in- 
ternal error estimate, e, while an external error estimate is 
provided by the standard deviation, tr, of repeated obser- 
vations at different epochs. From these and the number of 
observations, n, the probability P{x 2 ) that the observed 
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Fig. 3. Distribution of the mean errors of the mean radial 
velocities in the catalogue {top) , and the time span covered 
for each star (bottom). 



scatter is due to measuring errors alone may be computed 
as described in greater detail by Andersen & Nordstrom 
(1983b). P( X 2 ) < 0.01 is adopted as our criterion for cer- 
tain velocity variability - mostly due to binary orbital 
motion - for both the CORAVEL and CfA data. 

Normally, the mean error of the mean radial velocity is 
computed as a * n^ 1 / 2 . However, if fortuitous good agree- 
ment between a few observations results in a < e, then 
e * n -1 / 2 is given instead as a more realistic estimate of 
the mean error of the average velocity. 
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Occasionally, a double correlation peak may identify 
a spectroscopic binary from just a single observation, but 
normally two or more observations are available. In such 
cases, the centre-of-mass velocity and the mass ratio of a 
double-lined binary may be computed by the method of 
Wilson (1941) without a full orbital solution, if the ve- 
locities can be properly assigned to the two components. 
For the 510 systems for which this has been possible, the 
systemic velocity is given instead of the raw average of the 
observations, and the mass ratio is given in Table 2 (elec- 
tronic form only). The binary population of the sample is 
further discussed in Sect. 15.11 

Fig. |3| shows the distribution of the mean errors of the 
mean radial velocities in the sample (upper panel), and 
of the time span covered by the observations of each star 
(lower panel). As will be seen, the mean error of a mean 
radial velocity is typically ^0.25 km s _1 and only rarely 
exceeds 1 km s _1 . The observations typically cover a time 
span of 1-3 years, but occasionally extend over more than 
a decade. 
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3.3. Rotational velocities 

For stars with significant rotation, the width of the cross- 
correlation profile is a good indicator of vsini . For the stars 
with CORAVEL observations, vsmi has been computed 
using the calibrations of Benz & Mayor <|1980llT5%4")> . For 
the CfA observations, the i>sim of the best-fitting tem- 
plate spectrum is a good measure of the rotation of the 
programme star (Nordstrom et al. 119941 11997b(l . 

For the slowest rotators, more elaborate procedures are 
needed to derive very accurate rotational velocities; for the 
fastest rotators, the shallow cross-correlation profiles yield 
results of low accuracy. Accordingly, vsim as derived from 
the observations are only given to the nearest km s _1 , 
and from 30 km s _1 and upwards only to the nearest 5 
or 10 km s _1 . As seen in Fig. 01 the great majority of the 
programme stars have rotations below 20 km s _1 . 

3.4. Parallaxes 

Good distances are crucial in order to compute accurate 
absolute magnitudes, space motions, and parameters de- 
rived from them. Trigonometric parallaxes, generally of 
very good accuracy, are available from Hipparcos for the 
majority of our relatively nearby programme stars (ESA 
I1997J) . Fig. shows the distribution of the parallaxes (it) 
and their relative errors (ov/V); most are better than 10%, 
nearly all better than 20%. The computation of distances 
for all our stars is discussed in Sect. 14.41 

3.5. Proper motions 

Accurate proper motions are available for the vast ma- 
jority of the stars from the Tycho-2 catalogue (H0g et 
al. I2000|l . This catalogue is constructed by combining the 
Tycho star-mapper measurements of the Hipparcos satel- 



Fig. 4. Distribution of rotational velocities in the sample. 

lite with the Astrographic Catalogue based on measure- 
ments in the Carte du Ciel and other ground-based cat- 
alogues. By this procedure the baseline for determining 
proper motions is extended up to nearly a century, against 
only 3.5 years for the Hipparcos mission itself. 

For a few stars, mostly very bright stars or close bi- 
naries without a Tycho-2 proper motion, a Hipparcos or 
Tycho measurement has been used instead. The typical 
mean error in the total proper motion vector is 1.8 mil- 
liarcsec/year, corresponding to a mean propagated error in 
the space velocity of 0.7 km s _1 from the proper motions 
alone (i.e., neglecting parallax and radial- velocity errors). 

4. Derived astrophysical parameters 

In order for the stellar data to be useful in discussing 
the evolutionary history of the Solar neighbourhood, a 
number of astrophysically interesting parameters must be 
derived from the raw observational data. In most cases, 
calibrations of the photometric indices in terms of intrinsic 
parameters arc found in the literature, except as noted 
below. We discuss each calibration in turn in the following. 

4.1. Interstellar reddening 

E(b-y) can be computed for F stars with /3 observations 
from the intrinsic colour calibration by Olsen l|1988[l . It 
has been applied in the photometric temperature and dis- 
tance determinations if E(b-y) >0.02 and the distance is 
above 40 pc; otherwise the stars are assumed to be unred- 
dened. Most stars with no value of E(b-y) are late- type 
dwarfs within 40 pc, which will have negligible reddening 
anyway. As seen in Fig.EI very few of the stars have E(b- 
y) > 0.05 mag. We note that E(b-y) may be overestimated 
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for the hottest, brightest, and most distant early F stars 
CBurstein l2UB5|l . 



4.2. Effective temperatures 

Effective temperatures for all the programme stars have 
been determined from the reddening-corrected b — y, c%, 
and mi indices and the calibration of Alonso et al. I|1996[) 
which is based on the infrared flux method. The resulting 
temperatures have been compared to the determinations 



by Barklem et al. ( 2002), based on a fit to the Balmer line 
wings using the latest broadening theory. The agreement 



Fig. 6. Distribution of reddening values in the sample. 



is excellent, with a mean difference of only 3 K and a 
dispersion of 94 K. We have also compared our results to 
the spectroscopic excitation temperatures determined by 
Bensby et al. (2003) for 63 of our stars; the latter are on 
average 93 K higher than ours, with a dispersion around 
the mean of only 57 K. The distribution of the effective 
temperatures in the sample is shown in Fig. 
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Fig. 5. Distribution of Hipparcos parallaxes {top) and -5 
their relative errors (bottom) for the whole sample. 
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Fig. 7. Distribution of the sample in T e ff. 
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4.3. Metal abundances 

The accurate determination of metallicities for F and G 
stars is one of the strengths of the Stromgren uvby[3 sys- 
tem. Among the available calibrations, we have used that 
by Schuster & Nissen Q1989[l for the majority of the stars. 
In our sample, ~600 stars are covered by both the F and 
G star calibrations of Schuster & Nissen (198.9), and the 
mean difference in [Fe/H] is 0.06 with a dispersion around 
the mean of only 0.07. We have further compared these 
photometric metallicities with the homogeneous spectro- 
scopic values for F and G stars by Edvardsson et al. ( 1993 ) 
and Chen et al. J2000J). The agreement is excellent, with 
mean differences of only 0.02 and 0.00 dex and dispersions 
around the mean of 0.08 and 0.11, respectively. A further 
comparison with the compilation by Taylor H2003|) shows 
a mean difference of only 0.01 dex, but a larger disper- 
sion of 0.12 dex, as expected for a compilation from many 
sources of varying quality. 

Within the range of validity of the Schuster & Nissen 
( 1989) calibration, we thus find the photometric metallici- 
ties to have no significant zero-point offset and remarkably 
small dispersion when compared to high-quality spectro- 
scopic values. However, as pointed out most recently by 
Twarog et al. QZGGty . the Schuster & Nissen 1(1555) 1 cali- 
bration seems to give substantial systematic errors in the 
metallicity computed for the very reddest G and K dwarfs 
(6 — y > 0.46), where very few spectroscopic calibrators 
were available at that time. Because our sample contains 
an appreciable number of such red stars and more spec- 
troscopic metallicities in this range have become available, 
we decided to derive an improved metallicity calibration 
for these stars, as follows: 

From the high-resolution spectroscopic studies of 
Flynn & Morell l(T557)> . Tomkin & Lambert l(T555|l . 
Thoren & Feltzing (gDJDJ , and Santos et al. PH| . we 
have extracted metallicities for 72 dwarf stars in the 
colour range 0.44 < b — y < 0.59 and performed a new fit 
of the uvby indices to these values, using the same terms 
as the Schuster & Nissen l!1989fl G-star calibration. The 
resulting calibration equation is: 

[Fe/H] = -2.06 + 24.56mi - 31.61m? - 53.64mi(6 - y) 
+ 73.50m? (6 - y) + [26.34mi - 0.46ci - 17.76m?]ci 

The fit of the photometric metallicities from this cal- 
ibration to the spectroscopic reference values is shown in 
Fig. |SJ (open circles). The dispersion around the (zero) 
mean is 0.12 dex. 

Spectroscopic abundances for such cool dwarfs remain 
affected by both observational and theoretical uncertain- 
ties (see, e.g. Thoren & Feltzing 2000). However, the de- 
terminations selected here seem to represent the current 
state of the art, and we have used our new calibration to 
compute photometric metallicities for the ~1500 stars in 
our sample with b — y > 0.46. For the ~600 stars in the 
interval 0.44 < b — y < 0.46, the new calibration agrees 
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Fig. 8. Comparison between our final photometric metal- 
licities ([Me/H]) and the spectroscopic ([Fe/H]) values 
used to establish the calibrations. Open circles denote the 
cool (GK) stars, dots the hot (F) stars (see text). 



with that by the Schuster & Nissen l(T5%5)l to within 0.00 
dex in the mean, with a dispersion of 0.12 dex. 

About 2400 of our stars with high temperatures 
and low gravities are outside the range covered by the 
Schuster & Nissen l(T589) calibration. For these stars we 
have adopted the calibration of (3 and mi by Edvardsson 
et al. l|1993fl . when valid. For the stars in common, the 
two calibrations agree very well (mean difference of 0.00 
dex, dispersion only 0.05). For stars outside the limits of 
both calibrations, we have derived a new relation, using 
the same terms as Schuster & Nissen (1989) for F stars. In 
addition to the above new spectroscopic sources, we used 
Burkhart & Coupry JT55TJ), Glaspey et al. <(T553)) and 
Taylor (2003) to extend the coverage in b-y, mi, Ci, and 
[Fe/H]. From 342 stars in the ranges: 0.18 <b-y < 0.38, 
0.07 < mi < 0.26, 0.21 < ci < 0.86 and -1.5 < [Fe/H] < 
0.8, we derive the following calibration equation: 

[Fe/H] = 9.60 - 61.16mi + 81.25rm(6 - y) 
+ 224.65m?(6 - y) - 153.18mi(6 - y) 2 
+ [12.23 - 90.23mi + 38.70(6 - y)] \og(m x - c 3 ), 

where c 3 = 0.45 - 3.98(6 - y) + 5.08(6 - yf . 

The fit of these photometric metallicities to the spec- 
troscopic values is shown in Fig. |SJ (dots); the dispersion 
around the relation is 0.10 dex. For the stars in common, 
the new calibration and that by Schuster & Nissen l|1989f) 
again agree very well (mean difference 0.02 dex, dispersion 
only 0.04). More detail on the new calibration is given by 
Holmberg pTO)l . 
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Fig. 9. Distribution of metallicities for the whole sample 
(full histogram). For comparison, the dotted curve shows 
a Gaussian distribution with mean of -0.14 and dispersion 
of 0.19 dex, covering the same area as the histogram. 

The distribution of the photometric metallicities de- 
rived as described above is shown in Fig. [5] A Gaussian 
curve (with a mean of -0.14 and a dispersion of 0.19 dex) 
has been plotted to highlight the tail of metal-poor stars 
in the real distribution. This metallicity distribution for 
F- and G-type dwarfs is almost identical to the one found 
for K-type giants by Girardi & Salaris (2001), with a mean 
of -0.12 and a dispersion of 0.18 dex. 

4.4. Distances and absolute magnitudes 

Most of our programme stars are nearby and have trigono- 
metric parallaxes of excellent quality from Hipparcos (see 
Sect. 13.41 and Fig. We have therefore chosen to first 
determine distances for our stars based on the Hipparcos 
parallaxes, either directly or indirectly. The distances are 
used to compute tangential space motion components 
from the proper motions, and absolute magnitudes used 
in the determination of ages and masses. 

When the Hipparcos parallax is either unavailable or 
less accurate, a photometric parallax is used. We have 
adopted the distance calibrations for F and G dwarfs by 
Crawford 1)1975(1 and Olsen ©84); if both are valid for 
the same star, the F star calibration is preferred (Note 
that this calibration requires a /3 value) . We have checked 
the photometric distances against the subset of Hipparcos 
parallaxes with relative errors below 3% (Fig. HOfl . The 
trigonometric and (distance independent) photometric 
parallaxes agree very well, with no significant colour- 
dependent bias: The photometric distances have an un- 
certainty of only 13%. 
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Fig. 10. Photometric vs. Hipparcos distances for the sin- 
gle main sequence stars with parallax errors below 3%. 
Top: F dwarfs; bottom: G dwarfs. 

Accordingly, the Hipparcos distance is adopted if the 
parallax is accurate to 13% or better; otherwise we adopt 
the photometric distance. However, the photometric dis- 
tance calibrations are not valid for binaries, giants, and 
many kinds of peculiar stars. Such stars reveal them- 
selves by large discrepancies between the trigonometric 
and photometric distance estimates. The (few) stars with 
photometric distances deviating more than 3<r from the 
Hipparcos distances are nagged in the catalogue as sus- 
pected binaries or giants, and no photometric distance is 
given if the Hipparcos parallax is too inaccurate as a dis- 
tance indicator on its own (235 stars). Similarly, no dis- 
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tance is given for stars which lack the necessary photome- 
try (typically the /? index) and/or reliable Hipparcos par- 
allaxes or fall outside the photometric calibrations (1214 
stars). These stars are all quite distant and of marginal 
relevance to the overall sample. 

From the adopted distances and the observed V magni- 
tudes we have computed the absolute magnitudes given in 
the catalogue, correcting for interstellar extinction when 
known. Moreover, a SMy index has been calculated as the 
magnitude difference between the star and the theoretical 
ZAMS at the same colour and metallicity, as an indicator 
of the degree of evolution of each star. Fig. 1111 shows the 
distribution of My and SMy values for the sample. 




3.65 



3.85 3. 8 3.75 3.7 3.65 

Fig. 11. My (top) and 5My (bottom) vs. logT e // for our 
sample, which by design consists of F and G dwarf stars. 



4.5. Ages 

Individual stellar ages are crucial in order to place the 
observed chemical and kinematical properties of the stars 
in an evolutionary context. Because of their importance, 
we have devoted a great deal of effort to finding the most 
reliable way to determine ages for the stars in our sample 
and assessing their individual and systematic errors. 

We start by noting that observable diagnostics of stel- 
lar ages are basically: (i): chromospheric activity, and (ii): 
evolution away from the ZAMS in the HR diagram. Both 
have strengths and limitations, as discussed recently by 
Lachaume et al. (|1999|) and Feltzing et al. I|2()()1J> . 

Chromospheric age determinations rely on the decline 
of stellar activity with time (see Soderblom et al. [1991 
for an overview). The strength of this technique is that it 
can be used for both F, G, and K-type dwarfs, including 
very young stars. A drawback is that the chromospheric 
activity indicators (e.g. X-ray and Ca II emission) de- 
cay into invisibility at about the age of the Sun, so the 
method cannot be used for the older stars which are of 
main interest for Galactic evolution. A more basic prob- 
lem is the time variability of stellar activity, similar to the 
activity cycles and more dramatic phenomena of the Sun, 
such as the Maunder minimum. Further, stellar activity is 
caused mainly by rotation, which decreases with age but 
can be influenced by, e.g. tidal interaction in binary sys- 
tems (Kawaler 1989). A chromospheric age can therefore 
be completely wrong for reasons that cannot be clarified 
without additional (substantial) observational data. 

Isochrone ages are determined by placing the stars in 
the theoretical HR diagram (Fig. Ill|) , using the observed 
T e ff, M v , and [Fe/H] and reading off the age (and mass) 
of the stars by interpolation between theoretically com- 
puted isochrones. Edvardsson et al. I|1998J) exemplify this 
technique in the present context. Given the presence of ob- 
servational errors, isochrone ages can only be determined 
for stars that have evolved significantly away from the 
ZAMS, where all the isochrones converge. This precludes 
the determination of reliable isochrone ages for unevolved 
(i.e. relatively young) stars, and also for G and K dwarfs 
which evolve along the ZAMS for the first long period of 
their life. 

Many of our stars are considerably older than the Sun. 
Moreover, chromospheric activity indicators exist only for 
a small fraction of them. Accordingly, we have chosen to 
derive isochrone ages for our stars, recognising that mean- 
ingful results will not be possible for all our stars. 

4.5.1. Selection of stellar models 

Selecting an appropriate set of theoretical evolution mod- 
els and verifying its correspondence with the observed 
stars is the first crucial step in any determination of 
isochrone ages. The youngest stars in our sample arc mas- 
sive enough that the stellar models must incorporate con- 
vective core overshooting where appropriate. Several such 
models exist and are, in fact, very similar in the theoretical 
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plane, but employ rather different transformations to the 
standard colour systems (see, e.g., the detailed comparison 
in Nordstrom et al. I1997a|) . We have preferred, therefore, 
to compute effective temperatures and luminosities for the 
programme stars and compare with the models directly in 
the logT e // — M v plane. 

In preparation, we have compared the latest models 
from both the Geneva (Mowlavi et al. 119981 Lejeune & 
Schacrcr 2001) and Padova groups (Girardi et al. [2000 
Salasnich et al. 2000). The two sets of models yield essen- 
tially the same ages (to within 10%), but have significant 
and different limitations for our purposes: The Geneva 
models extend to very large ages, but are computed for a 
relatively coarse grid of masses > 0.8M Q , which precludes 
a proper determination of masses and ages for our coolest 
dwarf stars and leads to some numerical problems in the 
detailed isochrone interpolations. The Padova models ex- 
tend to stars well below the lower mass limit of our sample, 
but the isochrones are terminated at an age of 17.8 Gyr, 
complicating the proper computation of mean ages and 
age errors for the oldest stars in our sample. Ages much 
in excess of 17.8 Gyr exceed all recent estimates of the age 
of the Universe, however, so on balance, we have chosen 
the Padova models for the final age determination. 

4.5.2. Choice of model compositions 

The next issue concerns the choice of chemical com- 
position for the models. It has long been known (e.g 
Edvardsson et al. 113331 Reddy et al- ESPH that disk stars 
with [Fe/H] < exhibit an average enhancement of the 
a-elements which rises approximately linearly to [a/Fe] 
~ +0.25 at [Fe/H] = -1 and remains constant at that or 
perhaps a slightly higher level in even more metal-poor 
stars. The total heavy-element content of metal-deficient 
disk stars is thus somewhat higher than the heavy-element 
content of the Sun scaled by the observed [Fc/H]. 

Moreover, recent work (e.g. Fuhrmann[1998 , Bcnsby et 
al. l2003[) has found that thick-disk stars appear somewhat 
more a-enhanced than thin-disk stars in the range — 1 < 
[Fe/H] < 0, which spans the vast majority of our sample. 
There is, however, no consensus on a precise criterion to 
distinguish between stars of the thin and thick disks, in 
particular whether thick-disk stars are all extremely old 
and/or all moderately metal-poor. This makes it imprac- 
tical to identify the ~ 5% thick-disk stars and estimate 
separate a-enhancements for thin- and thick-disk stars. 
Moreover, while Padova isochrones are available for the 
Solar mixture of heavy elements as well as with an en- 
hanced a-element content for some values of [Fe/H], the 
assumed a-enhancement ([a/Fe] ~ +0.35) is considerably 
greater than appropriate for most of our stars. 

Fortunately, a simpler procedure appears sufficient. 
As demonstrated most recently by VandenBerg (2000), 
isochrones computed with solar-scaled and a-enhanced 
compositions are almost indistinguishable, provided the 
total heavy-element content Z remains constant. For the 
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metal-poor stars, we therefore select scaled-Solar compo- 
sition Padova isochrones with a somewhat higher [Fe/H] 
than the directly observed value, as described below. 

Specifically, we assume an a-enhancement that is zero 
for Fe/H] > 0, rises linearly with decreasing [Fe/H] to 
+0.25 dex at [Fe/H] = -1.0 and +0.4 dex at [Fe/H] = 
-1.6, then remains flat at +0.4 dex at all lower metal- 
licities. Following VandenBerg H2()()()|l . we then increase 
the observed [Fe/H] by 75% of the corresponding value 
of [a/Fe] as the best estimate of the total heavy-element 
content of each star. We note that VandenBerg (2000) 
found his procedure to be less satisfactory for the most 
metal-rich compositions, but those results were derived 
for a constant [a/Fe] = +0.3 dex, even for [Fe/H] > 0. 
The above approximation should remain fully satisfactory 
in the slightly metal-poor regime with the far smaller a- 
enhancements adopted here. 
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Fig. 12. Comparison between the observed stars (known 
binaries excluded) and Padova isochrones for 0, 2, 4, 6, 
8, 10, and 15 Gyr at [Fe/H] = 0.00 ± 0.02, -0.50 ± 0.05, 
-0.90 ± 0.20, and -1.50 ± 0.25, after allowing for the a- 
enhancement and temperature corrections discussed in the 
text. 



4.5.3. Adjusting the temperature scales 

Having described our procedure for choosing models of 
appropriate heavy-element content for stars of different 
metallicity, we ask whether these models provide a satis- 
factory fit to the observed stars. This is especially impor- 
tant for the effective temperatures, which are notoriously 
difficult to predict in an absolute sense from stellar mod- 
els as well as from observation (see, e.g. Lebreton l2001fl . A 
temperature mismatch between the models and observed 
stars will enter directly into the derived ages. 

Because our sample is expected to include very old 
stars, the comparison must be made on the unevolved 
main sequence, i.e. for M v > +5.5. Good agreement is 
found for Solar and very metal-poor compositions, but for 
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intermediate values of [Fe/H] the models are too hot by 
small, but significant amounts, as also found by Lebreton 
( 2001 ). Ignoring this offset would drive the low-mass stars 
to spuriously high ages, given the tight spacing of the 
isochrones in this mass range. 

Accordingly, we have applied temperature corrections 
to the models that amount to a <51ogT e // of -0.015 at 
[Fe/H] = -1.5, rising linearly to S\ogT eff = -0.022 at 
[Fe/H] = -1.0 and dropping linearly again to zero at [Fe/H] 
= -0.3. With these corrections, we obtain the isochrone 
fits to the lower main sequence shown in Fig. ^1 which 
we consider satisfactory. 

4.5.4. Statistical biases in age determinations 

The classical way to determine an isochrone age is to plot 
the observed stars and computed isochrones together in 
the theoretical HR diagram, either the log T e / / — My di- 
agram (Fig.ni top) or the logT e // vs. 5 My variety (Fig. 
im bottom) used by Edvardsson et al. 1)19931 see their 
Figs. 10-11). Errors are then derived by varying each of the 
independent variables logT e //, My, and [Fe/H] by their 
estimated observational errors and noting the changes in 
the resulting age. 

However, age is a highly non-linear function of posi- 
tion in the HR diagram; the distribution of the observed 
parameters is highly non-uniform as well; and the obser- 
vational errors are not always negligible compared to the 
ranges over which these distributions and the derived ages 
vary considerably. An age probability distribution func- 
tion computed without regard to these effects will there- 
fore be biased; moreover, in the simple, "classical" ap- 
proach it is also incompletely sampled. Biased ages and 
misleading error estimates are the likely result. 

Statistical biases effecting the determination of 
isochrone ages include the following: 

1. Stellar evolution accelerates strongly when stars leave 
the main sequence; therefore, the density of stars in the 
HR diagram will be much higher in the main-sequence 
region than away from it. This, in turn, causes more 
stars to be scattered by observational errors from the 
main sequence into the subgiant region than the re- 
verse, and leads to a bias in favour of high ages. Note 
that this effect is in fact exacerbated if stars with 
poorly determined ages are eliminated, since unevolved 
stars necessarily have poorly-determined ages. 

2. Standard initial mass functions (IMF) rise towards 
lower stellar masses; a given isochrone will therefore 
not have equal numbers of stars in equal mass steps, 
but the density of stars will rise towards the ZAMS. 
Ignoring this effect will also lead to a positive age bias. 

3. Standard disk metallicity distributions (see Sect. Ifi.ljl 
contain many more metal-rich than metal-poor stars; 
observational errors will therefore again scatter more 
stars from the metal-rich peak of the distribution into 
the metal-poor tail than in the opposite direction. This 
will cause the corresponding ages to be derived from 



isochrones that are too metal-poor, i.e. too hot, and 
again a positive age bias results (the converse argu- 
ment applies to the tail of "super metal- rich" stars). 
4. Apart from such "intrinsic" effects in the data, the 
distribution of stars in the HR diagram will be non- 
uniform because, e.g. of a non-uniform age distribution 
of the stars themselves, or as a result of the criteria 
used to define the sample. Notably, the distribution in 
the HR diagram of our full magnitude-limited sample 
will be quite different from that of the volume-limited 
subsample, due to the inclusion of luminous, evolved 
stars from much larger distances. 

4.5.5. Age determination for the sample stars 

The techniques we have developed to allow for these bi- 
ases are superficially similar to those discussed recently 
by Lachaume et al. (1999) and Reddy et al. 12003) as re- 
gards the treatment of the evolution bias referred to above. 
There are, however, important differences, in that we treat 
all three parameters logT e //, My, and [Fe/H] equally, in- 
clude several additional sources of bias, and consider the 
whole chain of astrophysical links from data to age. Our 
method is outlined below and described in greater detail 
by J0rgensen & Lindegren ( 2004 ) . 

Briefly, for every point in a dense grid of interpo- 
lated Padova isochrones we compute the probability P 
that the star could in reality be located there (and thus 
have the corresponding age), given its nominal position 
in the three-dimensional HR "cube" defined by \ogT e ff, 
My, and [Fe/H]. To do so, we assume that the associated 
observational errors have a Gaussian distribution: 

P = exp(-(AT e// ) 2 /24 e// )*exp(-(AM,) 2 /2a2 f J* 

ex P (-(A[Fe/if]) 2 /2af Fe/ff] ) 
Here, AT e // etc. are the differences between the observed 
parameters of the actual star and the isochrone points 
considered. We assume constant errors (<r) of 0.01 dex in 
logT e // and 0.1 dex in [Fe/H] throughout; for My we use 
the individual error estimate if an Hipparcos parallax bet- 
ter than 13% exists; otherwise, the standard photometric 
value of 0.28 mag in (m-M) is adopted. 

Integrating over all points gives the global likelihood 
distribution for the possible ages of the star, conditioned 
to account for observational biases as described below. 
We call this the "G-function" and normalise it to unity 
at maximum. The most probable age for the star is then 
determined as the value for which the G-function has its 
maximum (see Fig. 11311 . 

The determination of the maximum value itself is a 
non-trivial task. Because of numerical noise due to the 
finite sampling of the isochrones, a simple maximum of the 
raw function results in spurious high-frequency features in 
the derived age distributions which go undetected in small 
samples, but have dramatic effects in densely-populated 
diagrams such as Figs. 1271 and 1301 The median of, say, the 
upper 50% of the function yields a more stable estimate, 
but if the corresponding age range includes one of the 
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Fig. 13. Top: Three slices of the three-dimensional HR 
"cube" showing the observed point (central panel) and 
three cuts through the 1-cr error ellipsoid. The probabil- 
ity distribution function (G-function, bottom) is computed 
from all points on the isochrones, not just those on or in- 
side the error ellipsoid. 

limits (0 or 17.8 Gyr), the estimate will be biased away 
from the limit, leading to spuriously low ages for the oldest 
stars (cf. Fig. I14|) . A Gaussian fit as used by Reddy et al. 
( 2003) is also more stable, but is a poor approximation in 
the frequent cases when the G-function is distinctly non- 
Gaussian (Fig. [H|l. 



After extensive tests with simulated and real data, 
smoothing the G-function slightly with a kernel depending 
on the width of the unsmoothed function was found to be 
the optimum procedure. The maximum of the smoothed 
function then yields a stable age estimate without signif- 
icant bias. Fig. ^| illustrates the procedure in the well- 
behaved case of a star located in a region of the HR di- 
agram where the isochrones are well separated, and the 
maximum of the G-function yields a well-defined age. 

Finally, the maximum value of the probability function 
found for any point on the isochrones is a measure of the 
degree to which the star is covered by the models. A small 
value signifies peculiar stars or large observational errors 
which preclude any realistic age determination. 

In computing the G-functions, we have accounted for 
the biases described in Sect. I4.5l4l as follows: 

1. Speed of evolution. The isochrones are more widely sep- 
arated in phases of rapid evolution, so such phases au- 
tomatically receive lower weight in the integrations. 
We emphasize that not only points within a 1-cr (or 
3-cr, Reddy et al. I2003[) error ellipse are included, but 
all points on the isochrones in all three dimensions. 

2. Stellar mass function. The varying density of stars on 
an isochrone towards the main sequence is accounted 
for by weighting each point according to the IMF 
(Kroupa et al. 1993). It can be argued that the slope of 
the actual distribution in the magnitude-limited sam- 
ple will be lower than that of the IMF due to the pref- 
erential inclusion of brighter, higher-mass stars, but 
the effect is hard to quantify and in any case small, as 
the range in masses covered by the sample is small. 

3. Metallicity bias. The excess of apparently metal-poor 
stars caused by observational scatter from the large 
peak of stars of near-Solar metallicity is straightfor- 
ward in concept. Allowing rigorously for it in prac- 
tice is another matter: A fully Bayesian approach re- 
quires an estimate of the a priori distribution which 
is a priori unknown and, moreover, rather different 
for the complete, magnitude-limited catalogue and for 
the volume-limited sample which will no doubt be pre- 
ferred in many applications (compare Figs. l9l and 126(1 : 
other subsamples would no doubt be different again. 
It appears unreasonable that the catalogued age of a 
given star should depend on the subsample of stars 
discussed together with it. 

Moreover, the first-order astrophysical effects consid- 
ered earlier in the procedure already seem to allow for 
these effects. First, if significant, the metallicity bias 
should appear as an excess of positive residuals at low 
[Fe/H] when photometric metallicities are compared 
with spectroscopic determinations; this is not seen in 
our data (nor by Edvardsson et al. [1993). Second, 
our revised metallicity calibration (see Fig. [SJ , by de- 
sign, yields the correct mean spectroscopic [Fe/H] for 
a given mean photometric metallicity. 
Finally, and probably most importantly in view of the 
sensitivity of the derived ages to small temperature 
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shifts, the samples of stars used to "normalise" the 
temperature scale of the models to that of the observed 
stars (Fig. 112(1 will already be affected by any residual 
metallicity bias. Our temperature shifts will therefore 
allow for it to first order. In fact, it could be argued 
that these corrections may, if anything, be too large 
because the stars were drawn from the full, magnitude- 
limited sample which preferentially includes young, 
metal-rich stars. 

In summary we believe that, for general use, little if 
any metallicity bias of significance remains in the ages 
given in the catalogue. If particularly precise ages are 
needed for certain types of stars, well-defined subsam- 
ples should be extracted and all steps in the analy- 
sis reviewed and/or repeated, including the tempera- 
ture and metallicity calibrations, [a/Fe] ratio(s), model 
compositions and bolometric corrections, and the a pri- 
ori distributions of the relevant parameters. 
4. Age bias etc. A strongly peaked age distribution (e.g. 
due to a starburst) could give biases analogous to those 
discussed above. No such peak is expected, and its ef- 
fects would again depend on the (sub)sample consid- 
ered. We have also not imposed any upper limit on the 
derived ages: While true ages greater than ~13 Gyr are 
implausible, observational error will cause some deter- 
minations to exceed this limit, and imposing a cutoff 
will bias the mean age of the oldest stars. 

In cases (l)-(2) we correct for the a priori information 
in a fully Bayesian manner. Including also a priori metal- 
licity and age distributions in the age determination at 
this stage would build our prejudices concerning the en- 
richment and star-formation history of the disk into our 
age estimates for individual stars. We have therefore not 
made such corrections to the ages listed in the catalogue 
(equivalent to assuming flat a priori distributions.) 

Using the above procedures, G-functions have been 
computed for all stars in our sample with the necessary in- 
put data, including binary stars etc. An electronic table of 
these functions, which illustrate the determinacy of each 
age determination at a glance, will be made available to 
interested readers by request to the corresponding author 
(B.R.J. ; bjarne@astro.lu.se). 

Bayesian probability theory offers an independent, al- 
ternative method to compute unbiased age estimates un- 
der similar conditions, provided that the a priori distribu- 
tions of the relevant parameters can be estimated with suf- 
ficient accuracy. An end-to-end Bayesian approach of this 
type is described and applied to the data of Edvardsson 
et al. {MS} by Pont & Eyer piHlj) . 

4.5.6. Estimating errors for the ages 

Realistic error estimates are crucial in any applications 
of the ages. From a well-behaved G-function such as that 
shown in Fig. 1131 we derive (separate) 1-a lower and up- 
per age limits as the points where the G-function reaches 
a value of 0.6. Extensive Monte Carlo simulations using 



artificial stars with typical observational errors have con- 
firmed that indeed 68% of the recovered ages then fall 
within ±lcr of the correct age. 

When both the upper and lower 1-a age limits fall 
within the range of the isochrones, 0-17.8 Gyr, the cat- 
alogue lists both the most likely age of the star as well 
as its upper and lower limits. In the following, we refer 
to such cases as "well-defined" ages (note that this term 
by itself does not imply a small error, only that the error 
estimate is reliable!). For more demanding applications, 
the sample should no doubt be restricted to single stars 
and perhaps also to stars with age errors below a specified 
limit; all information needed to do so is readily available 
in the catalogue. 

If the G-function peaks within the valid age range but 
one of the limits is outside it, that limit is not given in the 
table, indicating that the age is uncertain and its error 
also poorly defined. For stars near or beyond the limits of 
the isochrone set (very young or very old stars with large 
observational errors, duplicity or other spectral peculiari- 
ties) , the G-function may peak at or even outside the age 
limits of the isochrones (see Fig. 114(1 . In such cases, no 
value is given for the age, only the estimated upper or 
lower limits. 




Fig. 14. Right: Examples of stars located in regions of the 
HR diagram where reliable ages cannot be determined, 
and (left) the corresponding G-functions. The stars plot- 
ted as plus, asterisk, triangle, and diamond symbols in the 
HR diagram correspond to the solid, dotted, dashed, and 
dot-dashed curves, respectively. Isochrones are for 0, 2, 4, 
6, 8, and 10 Gyr. 

Finally, for the lowest-mass stars which have not 
evolved perceptibly, the G-function will show no well- 
defined maximum (see Fig. 114(1 . If the G-function is too 
flat to reach the 1-a confidence level (0.6) anywhere in 
the range 0-17.8 Gyr, or if the maximum probability value 
entering the computation of the G-function indicates that 
the star falls significantly outside the isochrone set, no age 
is given at all. 
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Fig. 15. The distribution of upper (dotted line) and lower Fig. 16. Relative age error (mean of lower and upper 
(solid line) l-a relative errors for the ages in the catalogue, bounds) vs. age for the 11,445 stars with "well-defined" 



ages. 



The distribution of the relative age errors is shown 
in Fig. El while Fig. El shows the mean relative age er- 
ror as a function of age for the stars with "well-defined" 
ages. Note that the condition that both lower and upper 
age limits should be determined removes old stars from 
the upper right in Fig. 1161 The impression of increasing 
precision with age that results when a numerical cutoff 
is mistaken for a physical upper age limit (e.g. Fig. 4 of 
Feltzing et al. 200 1 i is, of course, an illusion. 

We have compared our error estimates with those de- 
rived in the classical manner, i.e. by varying each input 
parameter by ±ler and adding the age errors in quadra- 
ture. We find that the latter are often underestimated by 
almost a factor 2. We attribute this to three causes: (i): 
Varying only one parameter at a time significantly un- 
derestimates the true range of values over which the age 
variations must be explored; (ii): the technique effectively 
samples a total of only six points on the three-dimensional 
probability function which we integrate in detail to com- 
pute the G-function; and (Hi): the standard way to add 
errors assumes implicitly that the G-function is Gaussian, 
which is manifestly not the rule (cf. Fig.ll4[). Some of the 
error estimation techniques for isochrone ages in earlier 
literature have in reality estimated fitting errors rather 
than true uncertainties, and the errors have likely been 
significantly underestimated in several cases. 

In the total sample of 16,682 stars, these criteria yield 
age estimates for 13,636 stars (82%), of which 11,445 stars 
(84%) have "well-defined" ages by the above definition. 
9,428 (82%) of these well-defined ages have estimated er- 
rors below 50%, 5,472 (47%) even below 25%. Eliminating 
known binaries of all types leaves us with 9,158 presum- 
ably single stars (83% of all 11,060 such stars in the sam- 



ple) with derived age values, 7,566 (83%) of which have 
ages that qualify as "well-defined". Of these in turn, 6,144 
single stars (81%) have ages better than 50% and 3,528 
(46%) better than 25%, respectively. 

Fig. El shows the distributions of derived ages for 
the complete (magnitude-limited) sample for increasingly 
strict limits on the accuracy of the ages, and also compares 
the distributions of the magnitude-limited and volume- 
limited samples. We stress that, due to the biases in the 
selection and age computation procedures already dis- 
cussed, none of these diagrams has a simple interpreta- 
tion in terms of the star formation history of the Solar 
neighbourhood. 

4.5.7. Checking the results 

We have subjected our age derivation procedures to a wide 
range of numerical and other checks. In analogy with the 
Monte Carlo simulations used to verify the error estimates 
as described above, we have created artificial samples of 
stars of specified ages and a range of masses from the 
isochrones, computed uvbyfi indices and My values us- 
ing the reverse transformations of those applied to the 
observed data, and added realistic random errors to the 
simulated observations. These artificial stars have then 
been subjected to our age determination procedure in ex- 
actly the same manner as the observed stars. We find 
that, within the limits of the computed errors, we recover 
the input ages without significant systematic error (but of 
course with a loss of stars in the parts of the HR diagram 
where ages cannot be determined). 
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Fig. 17. Top: Age distributions for single stars in the 
full sample. The curves show (top to bottom): (i) all 
ages (solid); (ii) "well-defined" ages (dots); (in) ages with 
errors < 50% (short-dashed); and (iv) ages with errors 
< 25% (long-dashed). Bottom: "Well-defined" ages in the 
full, magnitude-limited sample (solid) and in the volume- 
limited subsample for d < 40 pc (dots). Note that increas- 
ing the demand on accuracy progressively removes most 
of the oldest stars. 



Another reality check is to determine ages for stars in 
open clusters with good uvby photometry and known red- 
dening (including, but not limited to such favourable cases 
as NGC 3680; Nordstrom et al. H997a|) as if they were sin- 
gle stars. Again, good agreement is found with the results 



of detailed isochrone fits to the entire cluster sequence, 
but of course the unevolved lower main-sequence stars 
yield large error bars. Similar consistency checks have been 
made in binary stars with good data for the individual 
components. Further discussion of these tests is given by 
J0rgensen & Lindegren 1|2004[1 . 

A particular concern was to ensure that our proce- 
dure does not introduce metallicity-dependent systematic 
errors that could distort the resulting age-metallicity rela- 
tions (AMR). In order to do so, we have created artificial 
AMRs of specified shape, with a distribution of metallici- 
ties at each age corresponding to the observational scatter, 
and with a uniform distribution of stellar masses from the 
ZAMS value to the maximum reached for the assigned 
age and metallicity. For the experiment, we assumed both 
an AMR with [Fe/H] increasing linearly with time and 
another (completely unphysical) AMR in which [Fe/H] 
decreased linearly with time. Simulated observations and 
realistic errors were computed and ages and metallicities 
rederived from the artificial data. As before, many stars 
were lost for which reliable ages could not be determined, 
but those with small calculated errors delineated the in- 
put AMR without any systematic error - a result which 
inspires confidence in our method. Details on these simu- 
lations are given in Holmberg ( 2004 ) . 

A further external check was made by comparing with 
the ages by Edvardsson et al. 1)1993[1 . Of their 182 stars 
with ages, our procedure yields estimates of any quality 
for 179 stars and "well-defined" ages (see above) for 160 
stars. Fig. ^| compares the two sets of ages for the latter 
sample; a linear fit yields the relation 
log Age Edv g 3 = (0.11 ± 0.03) + (0.89 ± 0.04) * log Age our . 
I.e., our ages are on average slightly smaller than those 
of Edvardsson et al. (1993 ) for younger stars (where their 
My derived from 8c\ is biased towards brighter values), 
while our ages agree well for the oldest stars. The scatter 
around the 45° line is 0.12 dex; the fit only reduces it to 
0.11 dex. Edvardsson et al. I|1993|l estimated that their 
ages had errors of about ±0.1 dex. Our estimated mean 
relative error for the same sample is also 0.10 dex, some- 
what larger for the youngest and smaller for the oldest 
stars, but the sample excludes the very oldest stars for 
which an upper age limit cannot be properly determined. 

The agreement between our ages and those derived 
by Edvardsson et al. (1993) is quite gratifying: One would 
have expected a larger dispersion in Fig. llSl if the errors in 
the two age determinations were completely uncorrelated, 
which suggests that the errors are primarily observational 
rather than systematic. We recall that even though the 
same photometric data have been used, the stellar models 
(including opacities and convection descriptions), temper- 
ature and metallicity calibrations, absolute magnitude de- 
terminations (Hipparcos parallaxes vs. 5My derived from 
the 8c\ index), treatment of the a-enhancement in metal- 
poor stars, and the method for computing the ages have 
all changed in the intervening decade. Finding a relation 
such as that seen in Fig. 1181 strengthens one's confidence 
in the whole procedure. 



B. Nordstrom et al.: The Geneva-Copenhagen survey of the Solar neighbourhood 



17 



20 



10 



6 - 



mm 



•••V 

9 m m 



• • • . 



4 6 8 

Age (Gyr) This work 



10 



20 



Fig. 18. Ages from Edvardsson et al. ( 1993 ) vs. our results 
for the 160 stars in common with "well-defined" ages. 

We have also compared our results with ages computed 
with the Bayesian method of Pont & Eyer l|2004[l ■ Apart 



from a scale difference of ~ 25%, due to a different metal- 
licity scale and their choice of the median rather than 
maximum of the probability distribution as the preferred 
age, there is no significant difference between, e.g., age- 
metallicity relations derived with the two sets of ages. 

For completeness, we finally compared our new age de- 
terminations (with errors <25%) with the chromospheric 
ages derived by Rocha-Pinto et al. f2T)DD|) for the 85 stars 
in common in the range 1-20 Gyr. The plot is quite simi- 
lar to Fig. 8 of Feltzing et al. 1|2001[1 . with no trace of any 
correlation between the two sets of ages. 

In conclusion, while noting the uncertainties, we con- 
sider our age scale to be the best presently available for 
the whole sample. However, if the age of a particular sub- 
group of stars, e.g. thick-disk stars, is needed to the high- 
est precision, then the detailed elemental composition of 
a carefully defined sample of such stars should be deter- 
mined and models of precisely the same composition be 
computed to derive their ages. 

It must be emphasised that our reliability tests per- 
tain to the ages derived for individual stars. For any re- 
alistic distribution of true ages for a complete sample of 
stars, one must take into account that many of the least- 
evolved stars, especially the old low-mass stars, will have 
no derived age in the catalogue, simply because the ob- 
servations cannot measure their evolution. We caution, 
therefore, that the true age distribution of the full sam- 
ple cannot be derived directly from the catalogue; careful 
simulation of the biases operating on the selection of the 
stars and their age determination will be needed to obtain 
meaningful results in investigations of this type. 



Finally, we recall that many stars in our sample are bi- 
nary or multiple systems, for which the derived ages (and 
metallicities) will be unreliable. In general, there is in- 
sufficient information available to recover the data for the 
individual binary components from the combined photom- 
etry, but the great majority of the visual and spectroscopic 
binaries in the sample are known and identified in the cat- 
alogue (see !5.1|) . and can thus be excluded in studies where 
absolute statistical completeness is not important. 



4.6. Masses 

Mass estimates are needed as the main clue to the evolu- 
tionary history of the stars in the sample. Because each 
point on a model isochronc corresponds to a specific mass 
value as well as an age, we can compute an M-function de- 
scribing the probability distribution of model masses for 
an observed star from the Padova models, exactly analo- 
gous to the G-functions for the ages (see Sect. W. 



The M-functions are much better behaved than the G- 
functions and generally yield good masses also for stars to 
which no meaningful age can be assigned. Individual er- 
ror estimates are also given for all masses in the catalogue; 
they average about 0.05 M©. Fig. HHI shows the distribu- 
tion of the derived masses in both the magnitude- limited 
and volume-limited samples. The low-mass limit at 0.65 
Mq reflects the red colour cutoff of our sample. 
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Fig. 19. The distribution of derived masses for single stars 
in the full (solid line) and volume-limited sample for d < 
40 pc (dotted line). 
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4.7. Space velocities 

Space velocity components (U, V, W) have been computed 
for all the stars from their distances, proper motions, 
and mean radial velocities. (U,V,W) are defined in a 
right-handed Galactic system with U pointing towards the 
Galactic centre, V in the direction of rotation, and W to- 
wards the north Galactic pole. No correction for the Solar 
motion has been made in the tabulated velocities. Our ra- 
dial velocities are of superior accuracy (Fig. |3J and the 
average error of the Tycho-2 proper motions corresponds 
to only 0.7 km s^ 1 in the tangential velocities, so the dom- 
inant source of error in the space motions is the distance. 
Accounting for all these sources, we find the average error 
of our space motions to be 1.5 km s _1 in each component 
{U, V, and W). 

Fig. HOI displays ( U, V, and W) for all stars in the sam- 
ple with a measured radial velocity. The [/-Wand V-W 
diagrams show a smooth distribution, the Hyades being 
the only clearly discernible structure (cf. Fig. ^ and Sect. 
12.2(1 . The U-V diagram, on the other hand, shows abun- 
dant structure in addition to the Hyades, with four curved 
or tilted bands of stars aligned along approximately con- 
stant ^-velocities. Note that the preponderance of stars 
from the two nearest spiral arms in our local sample leads 
to the well-defined limits of these structures seen both in 
the U-V diagram and in the V-W diagram (at V = -30 
and +20 km s _1 ). 

These structures, discussed e.g. by Dehnen (1998), 
Skuljan et al. l{T?99). have velocities resembling classic 
moving groups or stellar streams and have been named 
(top to bottom) : The Sirius-UMa, Coma Berenices (or lo- 
cal), Hyades-Pleiades, and £ Herculis branches. The loca- 
tion of the £ Herculis branch in the U, V plane raises the 
interesting prospect that kinematically selected samples of 
local thick-disk stars (e.g. Bensby et al. 2003) might con- 
tain an admixture of somewhat younger, more metal-rich 
thin-disk stars. 

Because these structures do not consist of coeval stars 
(see Fig. l30fl . they are not simply the remnants of broken- 
up systems such as classic moving groups, but could be 
produced by transient spiral arm structures (De Simone 
et al. I2004[) or result from kinematic focusing by non- 
axisymmetric structures of the Galaxy such as the bar 
(e.g. Fux 2(I(TTI and references therein). Our data will al- 
low more refined searches for detailed dynamical substruc- 
ture resulting, e.g., from past merger events (e.g., Chiba 
& Beers IMiOl Helmi et al. l2T)0"3f or from the dynamical 
effect of the bar. 



4.8. Galactic orbits 

From the present positions and space motion vectors of 
the stars, we can integrate their orbits back in time for 
several Galactic revolutions and estimate their average or- 
bital parameters. Correlating such Galactic orbits with the 
chemical characteristics and ages of well-defined groups of 
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Fig. 20. U-V, U-W, and V-W diagrams for all stars in the 
sample. 
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stars may yield insight of great interest into their origin; 
see e.g. Edvardsson et al. (|1993f) and many later papers. 

Before using the observed space motions for this task, 
they must be transformed to the local standard of rest. 
For this, we have adopted a Solar motion of (10.0, 5.2, 
7.2) km s" 1 (Dehnen & Binnev HM^i . For the orbit in- 
tegrations we used the potential of Flynn et al. (TX996), 
adopting a solar Galactocentric distance of 8 kpc. The 
key properties of this model are: Circular rotation speed 
220 km s _1 , disk surface density 52 M Q pc~ 2 , and disk 
volume density 0.10 Mq pc~ 3 , in agreement with recent 
observational values (Reid et al. 119991 Backer & Sramek 
Flynn & Fuchs lT5Ml Holmberg & Flvnn ETftftHl . 

In the catalogue we give the present radial and vertical 
positions (R ga i and z) of each star as well as the computed 
mean peri- and apogalactic orbital distances R m in and 
Rmax, the orbital eccentricity e, and the maximum dis- 
tance from the Galactic plane, z max . Fig. 1211 summarises 
the distribution of these parameters. 

The detailed shape of especially the Rmin distribution 
is strongly dependent on the exact value of the Solar V- 
velocity applied, and the peaks can be related to the en- 
hanced density of stars in the streams seen in Fig. l2T)l 
Sirius-UMa causes the peak at ~ 8kpc in the R m in distri- 
bution, Hyades-Pleiades the one at ~ 7kpc, and the small 
bump at ~ 5.5 kpc can be traced to the C Herculis stream. 
Tests have shown that an increase of the Solar V-velocity 
by a few km/s changes the double-peaked structure into 
a smooth increase, due to the changed orbit distribution 
across the circular velocity. 




Fig. 21. Key parameters of the Galactic orbits for the 
whole sample. 



5. Statistical properties of the sample 

Few if any complete and unbiased samples of stars 
with complete, homogeneous astrophysical data exist in 
Galactic astronomy. Indeed, the conclusions of an obser- 
vational study may well depend more on the way the stars 
are selected than on the actual data obtained. Thus, for 
any application of the data presented here it is crucial 
to be aware of the manner in which the stars have been 
selected and observed, and how their astrophysical param- 
eters have been derived. 

A chief concern at the start of the project was to avoid 
any kinematic bias in the selection of stars. Therefore, 
stars for the photometric surveys were originally selected 
based on the HD visual magnitudes. These are known 
to have appreciable errors, so some stars brighter than 
the limits described in Sect. 12. ll will certainly have been 
missed, while many of the newly measured magnitudes 
were fainter than those limits. The wide range of HD spec- 
tral types included in the surveys should ensure that few 
if any metal-poor stars were missed from the final (much 
smaller) sample of F and G dwarfs defined from the mea- 
sured uvby(3 indices. The cool limit is somewhat less rig- 
orously denned, as old metal-rich dwarf stars might be 
mistaken for giants because of their strong CN bands, pos- 
sibly even in the MK spectral classes used to extend the 
red limit of the sample to K2 V south of S = -26°. 

As a practical factor, many of the radial- velocity obser- 
vations were made from preliminary versions of the pho- 
tometric catalogues, so the final selection criteria did not 
necessarily coincide precisely with the way the observa- 
tions were made. Stars that had already been observed 
could be excluded in the final sample if needed, but new 
stars could not necessarily be added. This is the case for a 
fraction of the unevolved F stars; they had been omitted 
initially because no reliable ages can be determined for 
them, and not all were included in the Hipparcos Input 
Catalogue which was later used to complete the data. 
Also, several hundred stars yielded no radial-velocity de- 
termination because of fast rotation, bright hot compan- 
ions, or for other reasons. 

Inevitably, when many tens of thousands of observa- 
tions are performed manually on hundreds of nights, occa- 
sional misidentifications and other errors occur. The thou- 
sands of double stars of all possible varieties provide rich 
opportunities for ambiguity and inconsistency in the at- 
tribution of individual observations. Great effort has been 
devoted to detecting and eliminating errors of all conceiv- 
able kinds, but not all cases can be resolved and rejected 
observations could not always be repeated. 

Finally, we caution again that the astrophysical pa- 
rameters for any star in the catalogue flagged as a binary 
are likely to be inaccurate and potentially misleading. 

5.1. Binary stars 

Binary stars are abundant amongst field stars, especially 
in samples limited by apparent magnitude because binary 
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stars are on average brighter than single stars. If unrecog- 
nised, the binary stars will appear as a substantial, but 
unknown fraction of the sample for which the derived 
ages, metallicities, distances, and space motions will all 
be wrong. Therefore, while corrected values cannot always 
be derived, identifying such cases is important, and much 
effort has been spent to that end. 

Our sample contains numerous double stars with sep- 
arations ranging all the way from invisible to fully re- 
solved in both photometric and radial-velocity observa- 
tions. Similarly, brightness ratios range from unity to sev- 
eral magnitudes, with associated observational difficulties. 
As far as possible, all such cases have been recorded at 
the telescope (see comments in the photometry papers), 
and stars with visual companions are flagged in the cata- 
logue and the available information given. As noted above, 
many radial-velocity observations of such visual compan- 
ions have been made; these will be made available sepa- 
rately at the web site of Observatoire de Geneve. 

The accuracy of our new radial velocities varies with 
the rotational line broadening, but the large majority are 
of excellent quality (cf . Fig. |3J) . Our average of more than 
four observations per star should allow good binary de- 
tection statistics, although the actual number varies con- 
siderably from star to star and many stars have only two 
observations. Accordingly, some long-period and/or low- 
amplitude binary stars will likely remain in the sample, 
but a large fraction of the systems with periods below 1000 
days should have been revealed by our data (see below). 
Many double-lined systems also manifested themselves by 
a double cross-correlation peak already at the telescope. 
All available spectroscopic information on duplicity for the 
stars in our sample is also given in the catalogue. 

The complete sample of 16,682 stars contains 3,537 
(21%) visual double stars and 3,223 (19%) spectroscopic 
binaries of all kinds. The total number of binary stars of 
any type is 5,622 (34%), as some visual binary components 
are also spectroscopic binaries. This fraction corresponds 
well to the total frequency (32%) of binary systems of all 
types with periods less than 10 5 days found for G dwarfs 
by Duquennoy & Mayor Ijl991|l . after careful correction 
for detection incompleteness. Our sample thus comprises 
11,060 stars that are not known to be double; of these, 
7,817 have measured radial velocities consistent with their 
being true single stars. 

A priori, we would expect a larger fraction of bi- 
nary systems in our magnitude-limited sample than in the 
volume-limited sample of Duquennoy & Mayor fUSEpI be- 
cause binaries are, on average, brighter than single stars. 
However, restricting the statistics to the volume-limited 
subsample within 40 pc, we still find a total binary frac- 
tion of 32%, suggesting that the great majority of binaries 
in the sample have indeed been detected. 

Finally, we warn that the mass ratios for 510 double- 
lined binaries derived from our (generally sparse) radial- 
velocity observations and listed in Table 2 should not be 
used to derived mass ratio distributions or any other sta- 
tistical properties for binary stars in general. This sam- 



ple is heavily biased towards binary stars with near-equal 
components and periods of a few days, which show signifi- 
cant line separation without excessive rotation (i.e. larger 
than ~30 km s _1 ). 

5.2. Magnitude completeness 

The distribution of the photoelectric V magnitudes for 
all stars in the sample as well as for those with complete 
astrophysical data are shown in Fig. From this dia- 
gram, the magnitude-completeness of the sample can be 
estimated by comparing with the distribution expected for 
a uniform volume density of stars. The full sample begins 
to depart from completeness near V = 7.6. 

The original surveys were designed with variable mag- 
nitude limits in order to better approximate a volume- 
complete sample. The magnitude limit thus depends on 
colour, as shown in Fig. [23 Note that a substantial frac- 
tion of the earliest F dwarfs (0.200 <b-y < 0.300) lack 
radial-velocity data due to fast rotation, while the cooler 
stars are essentially complete in this regard. 

The magnitude completeness as a function of colour 
may be characterized by two numbers: Vu m , the magni- 
tude at which incompleteness sets in, and V cu t , the magni- 
tude beyond which a negligible fraction of the stars were 
measured. Estimates of these completeness limits are given 
in the table below. 

vlim Vcut 

0.205 < b~^j < 0.300 7/7 8~1T 

0.300 < b-y < 0.344 7.8 8.9 

0.344 < b-y < 0.420 7.8 9.3 

0.420 < b-y < 0.540 8.2 9.9 



We note again that the coolest dwarfs are included only 
for S < —26°. These 1,277 stars are flagged in Table 1, so 
the sample can be cleanly divided into one subsample that 
is homogeneous over the sky to a constant colour limit, and 
another that is similarly homogenous to a redder colour 
limit, but covers only the southernmost 28% of the sky. 

5.3. Volume completeness 

The degree of volume completeness of the sample can be 
estimated from the computed distances of the stars. The 
programme stars have a wide range of absolute magni- 
tudes, but are drawn from the apparent-magnitude lim- 
ited photometric surveys, so the distance to which the 
sample is volume complete will vary considerably through 
the sample, primarily - but by no means exclusively - with 
the b-y colour index. 

Fig. El shows the distribution in distance of the four 
colour intervals defined in Fig. |23 As expected, the earli- 
est and intrinsically brightest F stars arc complete to the 
largest distance, ~70 pc (but note that many of these have 
no measured radial velocities). The G5 dwarfs as well as 
the later-type stars south of 5 = -26° appear to be com- 
plete to ~40 pc, as intended. We recall that 1,449 stars 
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Fig. 22. Distribution of observed V magnitudes for the 
full sample (full line), and for the stars with measured 
radial velocities (dotted line). The dashed curve shows the 
expected relation for a uniform, volume complete sample. 




Distance (pc) Distance (pc) 

Fig. 24. Distance distributions by interval in b-y colour 
(full line). The dashed curves show the expected relation 
for a uniform, volume complete sample. In the last panel, 
the lower histogram shows the distribution of stars south 
of S = -26°. 
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Fig. 23. Distribution of V magnitudes by interval of b-y 
colour. Markings as in Fig. 1221 



have no distance listed in the catalogue because of inade- 
quate data (see Sect. 14. 4|) . 

In the end, the volume within the limiting distance of 
40 pc contains only 1,685 stars with complete data, out 
of the 16,682 stars in the full sample (or the over 30,000 



stars in the original uvby surveys!). However, as distances 
are known for nearly all the stars, the absolute-magnitude 
bias is quantifiable and can be allowed for in computing 
the relative frequencies of the different types of star. 

5.4. Completeness of masses 

Stellar mass varies much more slowly over the HR dia- 
gram than age. Accordingly, the M-functions used to de- 
rive masses and their errors for our stars (Sect. I4.6|l are 
much better behaved than the corresponding G-functions 
for the ages, and useful mass estimates can be derived 
also for stars to which no meaningful age can be assigned. 
Accordingly, masses have been derived for all 14,381 stars 
with determinations of logT e //, My, and [Fe/H] and lo- 
cated inside the region in the HR diagram covered by the 
theoretical isochrones (see Fig. I19JI . 

5.5. Completeness of ages 

In Sect. 14. 531 we discussed how we have determined errors 
for the ages derived for the individual stars. The degree 
of completeness of the age data depends, then, on how 
large errors one is willing to accept. For the total sample, 
publishable ages as defined in Sect. 14. 5^1 are available for 
13,636 stars or 82% of the full sample (9,158 or 83% of 
the single stars), of which 11,445 (and 7,566 single stars 
- 69% of the total in both cases) have "well-defined" ages 
by the definition in Sect. 14.5.61 If the relative error limits 
are set to 50% or 25%, the numbers are 9,428 and 5,472 
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Fig. 25. Age vs. mass for stars with well-defined ages. 
The apparent correlation of age and mass is an artifact 
caused by the over-representation of evolved stars and the 
fact that reliable ages cannot be determined for unevolved 
stars. 

(57% and 33%), respectively, for all stars, and 6,144 and 
3,528 (or 56% and 32%) for the single stars. 

The fraction of stars for which reliable ages can be de- 
rived varies strongly over the HR diagram. For the evolved 
F stars the situation is relatively favourable. For stars on 
or near the ZAMS, only an upper limit to the age can, 
in effect, be determined. And for the still-unevolved low- 
mass stars no usable ages can be determined at all; in our 
sample, essentially no star below 0.90 M Q or fainter than 
My — 4.5 yields a meaningful age. 

These effects are evident in Fig. |2H1 which shows the 
correlation between the ages and masses derived for the 
stars in the sample. Three effects are obvious from the 
figure: (i): The upper envelope of masses for a given age, 
set by the evolution of the stars onto the giant branch 
(but blurred by metallicity effects); (ii): the lack of stars 
younger than ~ 1 Gyr, a direct result of our blue colour 
cut-off at b-y = 0.205; and (in): the lack, at each age, 
of the low-mass, unevolved stars for which well-defined 
ages cannot be determined. The few scattered stars below 
the main relation are probably undetected giants or other 
peculiar stars. 

6. Discussion 

The data base presented here provides a basis for a thor- 
ough re-evaluation of some of the global properties of 
the Solar neighbourhood and the Galactic disk. Several 
of these analyses will require detailed simulations of the 
predictions of various competing models, subjecting the 



simulated stellar populations to the same selection crite- 
ria as used to establish our samples, and comparing the 
results with the real data set. 

Such simulations will be the subject of future papers, 
but are beyond the scope of the present discussion. Here, 
we will simply review a few of the 'classic' diagnostic re- 
lations for the Galactic disk, i.e, the metallicity distribu- 
tion for long-lived stars, the age-metallicity relation for 
the Solar neighbourhood, the radial metallicity gradient, 
and the age-velocity and metallicity-velocity relations. 

In each case, two features are essential: (i) the lack of 
kinematic selection bias in our sample; and (it) the new 
radial- velocity data which allow to identify stars that have 
not taken part in the evolution of the local disk, but just 
happen to pass through it at this time; remove contami- 
nation by unrecognised binary stars; complete the veloc- 
ity information for the sample, and estimate scale-height 
corrections. At the same time, the very strong absolute- 
magnitude bias in our sample and limitation of the coolest 
stars to the cap south of —26° declination must be kept 
firmly in mind and corrected for as appropriate. 

As a preliminary, we warn that the true uncertainty 
of the ages of especially the oldest stars must be kept in 
mind when interpreting the following diagrams (note in 
this connection that the final refinements in our age com- 
putation methods have resulted in significant changes in 
the ages of the oldest stars since our preliminary discus- 
sion m Holmberg et al. E^UU^ . 

As shown in Fig. 1171 restricting the sample to stars 
with increasingly better-determined ages preferentially 
eliminates the oldest stars and makes it increasingly dif- 
ficult to determine the maximum age of stars in the disk. 
Conclusions based on single points in, e.g. Figs. |2"TI or |3*U1 
will be risky at best. E.g., we note that the statement by 
Sandage et al. (|2003jl . that "the age of the field stars in 
the solar neighbourhood is found to be 7.9 ± 0.7 Gyr", is 
based on a few stars assumed (but not known) to be super- 
metal-rich and ignores the existence of binary stars on the 
main sequence. At the same time, the authors recognise 
that the open cluster NGC 6791, with [Fe/H] = +0.3-0.4, 
is about 10 Gyr old, a warning against quick conclusions 
drawn from a single diagram. 

6.1. Metallicity distribution (the "G dwarf problem" ) 

Relative to the predictions of a closed-box model with con- 
stant initial mass function (see e.g. Pagel lTMTjl . the ob- 
served metallicity distribution function for unevolved low- 
mass dwarfs shows a significant lack of low-mass metal- 
poor stars (van den Bergh IT962I Schmidt 1963). These 
should have been formed along with the high-mass stars 
that produced the heavy-element content of later stellar 
generations, but have not been found in the predicted 
numbers in previous surveys ("the G-dwarf problem"). 
While closed-box models are now obsolete, the metallic- 
ity distribution of long-lived stars remains a fundamental 
constraint on any of its successors. 
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The G-dwarf problem could in principle have two dif- 
ferent explanations: (i): The 'missing' metal-poor dwarfs 
do in fact exist, but have not been found by previous 
searches; a population of slow-moving metal-poor stars 
could exist which would be missed in proper-motion sur- 
veys, but be included in our sample. Or (ii): the Galactic 
disk did not evolve according to simple closed-box models. 

The G-dwarf problem was rediscussed exhaustively by 
J0rgensen ()2000[l . based on an early version of the present 
data set. We refer the interested reader to that paper for 
the full discussion, since no material changes have been 
made to the underlying data since that study. 

Here we just recall the salient conclusion, i.e. that the 
solution to the G-dwarf problem is not a previously undis- 
covered metal-poor population. On the contrary, when 
drawn from a volume-complete and kinematically unbi- 
ased sample, the true fraction of metal-poor dwarfs is in 
fact only about half as large as previously believed, and 
the agreement with closed-box models even worse than 
before. Better physical models, not better data, are now 
the most urgent requirement. Fig.l26lshows the metallicity 
distribution for the volume complete part of our survey. 
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Fig. 26. Distribution of metallicities for the volume com- 
plete sample of single stars (full histogram). For compari- 
son the dotted curve shows the reconstructed distribution 
for G dwarfs from J0rgensen ; 2000;. which is corrected for 
scale height effects and measurement errors. 



6.2. Age-metallicity relations 

Since the pioneering studies of Mayor <|1974|l and Twarog 
(1980), the relationship between average age and metallic- 



ity in the solar neighbourhood has been a subject of con- 
tinuing debate, summarised in the recent study by Feltzing 
et al. I)2UU1|I . The debate concerns both the overall shape 
of the mean relation, whether there is scatter in the rela- 
tion over and above that due to observational errors, and 
if so, what might be the cause of this scatter. 

The accurate spectroscopic metallicities of Edvardsson 
et al. Q1993[l established beyond reasonable doubt that 
real scatter exists in the relation, far above that caused 
by observational errors. Little if any variation of mean 
metallicity with age was found, except for the very old- 
est stars (> 10 Gyr) where a downturn was observed in 
agreement with previous studies. However, as pointed out 
by Edvardsson et al. (1993) themselves, the restriction of 
their sample to F-type dwarfs automatically excluded any 
old, metal-rich stars if such existed, a fact that has been 
overlooked in several later discussions. The very recent 
similar analysis by Reddy et al. l|2003fl of a larger number 
of elements in 181 F and G dwarfs reached substantially 
the same conclusions as Edvardsson et al. (1993 ). 
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Fig. 27. Age - metallicity diagram for 7,566 single stars 
with "well-defined" ages in the magnitude-limited sample. 
Note that individual age errors may still exceed 50% (cf. 

Fig. Eg. 

The selection criteria underlying the present sample 
were specifically designed to include the old, metal-rich 
stars that would have been missed by Edvardsson et al. 
(1993) even if they existed. At the same time, it was re- 
alised that these are just the stars for which the derived 
ages will be the most uncertain, as demonstrated implic- 
itly by Feltzing et al. (2001). In this context we note that 
Reddy et al. ( 2003 ) find again a more clear-cut rise of mean 
metallicity with age than the two studies just mentioned, 
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Fig. 28. Age - metallicity diagrams for single stars in the 
volume-limited subsample within 40 pc; large dots show 
mean ages and metallicities in 10 bins with equal numbers 
of stars. Top: The 462 stars with "well-defined" ages (no 
limit on actual errors), bottom: The 142 stars with ages 
better than 25%; average error bars (14%) in three age 
bins are shown. 



primarily because the Reddy et al. sample contains few 
old, metal rich stars. It is unclear whether this is caused 
by their selection procedure, and it is therefore of inter- 
est to re-examine the age-met allicity diagram constructed 
from the stars with well-defined ages in the unbiased sam- 
ple presented here (see Figs. E3 and EEJ)- 



The most obvious features of Fig. [Z7\ remain direct re- 
sults of our selection criteria: The absence of stars near 
the lower left edge of the diagram is caused by our blue 
colour cutoff. The predominance of young, metal-rich stars 
is due to their intrinsic brightness and the correspond- 
ingly large volume they sample; the few young, apparently 
super-metal-rich stars are probably at least in part distant 
giant stars for which the interstellar reddening has been 
overestimated (see below and Burstein 2003). 

Apart from these features, little variation in mean 
metallicity is seen, except possibly for the very oldest 
stars which in general have kinematics characteristic of 
the thick disk. Even some of these have solar-like metal- 
licities, however, and we note that Bensby et al. (2003) re- 
cently derived a metallicity distribution for the thick disk 
that extends to stars with [Fe/H] > and Sun- like [a/Fe] 
abundance ratios; those stars are, however, also relatively 
young, so their thick-disk pedigree may remain open to 
question. The "£ Herculis branch" of disk stars in the U- 
V diagram fDehnen 119981 Skuljan et al. t!999j> could be a 
source of such stars. 

When interpreting Fig. |2"7I fand Fig. l5U|l . the substan- 
tial errors of even the "well-defined" ages should always 
be kept in mind; the presence of stars appearing to be as 
old as 14 Gyr is easily explained by observational errors. 
Uncertainties in the temperature scales of the observed 
stars and theoretical isochrones (Sect. |4~5~3*)) remain a po- 
tential source of systematic error, and numerical details 
of the age computation may introduce spurious features 
in diagrams such as Figs. E3 an d USB which can appear 
dramatic without the elaborate precautions described in 



Sect. 14.5.51 but may remain in more subtle form. 

In order to avoid the strong absolute-magnitude bias in 
Fig.EZI we plot in Fig.l2"51the same data for the volume- 
limited subsample within 40 pc. Despite the drastic re- 
duction in number of stars, the lack of young metal-poor 
stars produced by our blue colour cutoff remains well vis- 
ible and is responsible for the upturn of the mean relation 
for the youngest ages. The disappearance of the young 
"super-metal-rich" stars supports their interpretation as 
an artifact of the de-reddening procedure for distant giant 
stars. The lack of any overall metallicity variation in the 
thin disk is even more pronounced than before. Finally, 
the scatter in [Fe/H] at all ages again greatly exceeds the 
observational error of ~0.1 dex. Limiting the sample to 
the best-determined ages leads to the same conclusions. 

This picture of the age-metallicity distribution for field 
stars agrees well with the most recent studies of open clus- 
ters by Friel et al. PUU2|1 and Chen et al. (|2TTU3|l . These 
studies show the same constant mean metallicity and large 
scatter at all ages for the clusters as our study does for 
the field stars. The discussion by Sandage et al. (2003) 
also focuses on the existence of old, metal-rich subgiants 
in their sample (the 10-Gyr-old metal-rich cluster NGC 
6791 being another example) - as in fact pointed out al- 

mm 



ready by Stromgren (1963,). Clearly, more realistic models 
of the true complexities of star formation and chemical 
enrichment of the interstellar medium are required. 
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We recall that we have taken special pains to verify 
that our age computation technique will not distort the 
overall trends in the resulting age-metallicity diagrams (cf. 
Sect. l4.5.7Jl . Individual ages shown in Figs. 1771 andl^Slftopl 
may still be uncertain by 50% or more, however (cf. Fig. 
I15|) . which must be taken into account in any discussion 
of these diagrams. 

"Cleaner" subsamples of stars can be selected from 
the catalogue (e.g. Fig.051 bottom), but introduce further 
strong selection effects; cf. the top and bottom panels of 
Figs. EHl and 1201 Evaluating the true interplay of cosmic 
scatter in the chemical evolution of the disk with observa- 
tional errors and uncertainties in the age determinations 
will require detailed numerical simulations, which must in- 
clude models of the statistical biases discussed above and 
by Pont & Eyer ( 2004 ) . Such detailed analyses are beyond 
the scope of the present paper. 

6.3. Radial metallicity gradients in the disk 

In addition to the evolution of metallicity with age, our 
sample can be used to study the radial metallicity gradient 
in the Galaxy. From the R m in and R m ax of the stellar or- 
bit, the mean orbital radius R m can be calculated. Fig. [2^1 
shows the radial metallicity gradient for three groups of 
stars of different ages. The slopes of the fitted lines are 
-0.076±0.014, -0.099±0.011, and +0.028±0.036 dex/kpc. 
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Fig. 29. Radial metallicity gradient for single stars in 
three age ranges. R m is the mean radius of the stellar 
orbits. 

This can be compared to other studies using FG dwarfs 
and giants by Mayor fTHZSl, Cepheids by Andrievsky et 
al. I)2002|l . planetary nebulae by Maciel et al. I|2003l) . and 



open clusters by Friel et al. (|2002|l and Chen et al. 
In general they find the radial gradient to evolve over time 
from values between -0.02 and -0.06 for the youngest stars 
to between -0.08 and -0.12 for the older stars. This is com- 
patible with our two younger age groups which show a 
mild steepening with age of the radial gradient. The old- 
est stars in our sample, on the other hand, show no radial 
gradient at all. This is an indication that these stars do 
not follow the general evolution of the younger stars in the 
disk, but are of a different origin, perhaps from the thick 
disk. 

The radial metallicity gradient of ~-0.09 dex/kpc seen 
for stars younger than 10 Gyr can be used in an attempt 
to correct the age-metallicity distribution for the effects of 
radial migration in the disk (cf. Wielen et al. 11996(1 . This 
correction has no discernible effect on the distribution, and 
the dispersion around the mean for the volume-limited 
sample is unchanged: Computing the mean and standard 
deviation of the metallicities in the age range 2-12 Gyr in 
Fig. Hi (top), we find <[Fe/H]> = -0.16, a = 0.20 dex 
with no correction, <[Fe/H]> = -0.19, a = 0.20 dex after 
applying the correction for the radial metallicity gradient. 

Stars migrating into the sample from orbits centered 
elsewhere in the disk obviously account for only a minute 
part of the scatter in the age-metallicity diagram, and 
both dispersions are clearly much larger than can be ex- 
plained by errors in the metallicities, as also concluded 
by Edvardsson et al. ijlMSjl and Reddy et al. (|2TJU5|l from 
detailed spectroscopy. Our removal of binaries from the 
sample and careful study of the age uncertainties also 
guarantee that poorly-determined ages do not affect this 
conclusion (and the main trends of Figs. I28| are insensitive 
to even substantial horizontal redistribution of the points 
in the range 2-12 Gyr). 

6.4. Age-velocity relations 

The observed space velocity components U, V, W for all 
single stars in the sample are shown as functions of age in 
Fig. |20| with weak as well as strong limits on the accuracy 
of the ages. The population of bright, early F-type stars 
is prominent, but otherwise the diagrams reflect the slow 
increase of the random velocities with age which is at- 
tributed to heating of the disk by massive objects such as 
spiral arms or giant molecular clouds. The rate of change 
and maximum velocity dispersion reached in the thin disk 
are of key importance for the interpretation in terms of 
the local dynamics of the disk itself. 

In order to quantify these relations and also search for 
a kinematic signature of the thick disk, we have computed 
the velocity dispersions of the single stars with the best- 
determined ages. Fig.l31lshows the resulting relations and 
the power laws fitted to them, excluding the youngest and 
oldest bins. The resulting exponents are 0.31, 0.34, 0.47, 
and 0.34 for U, V, W, and the total velocity dispersion, re- 
spectively, with an uncertainty of 0.05, in close agreement 
with Holmberg ijM)!"]) and Binney et al. (200QJ>. Aside 
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Fig. 30. U, V and W velocities as functions of age for 
all single stars with complete data. Top panels: All 7,237 
stars with "well-defined" ages and all velocity data; bottom 
panels: The 2,852 stars with ages better than 25%. 



from the low exponent for the total velocity dispersion, 
the evolution of the velocities over time is characterised 
by a small increase in the ratio of ay / ajj and a larger in- 
crease in the ratio aw / &u . The consistency of these recent 
studies, which all use isochrone ages but otherwise differ- 
ent samples and methods, is in sharp contrast to studies 
using ages derived from chromospheric emission, which 
give values of either 0.26±0.01 or 0.59±0.04 (Hanninen & 
FlynnHHUl. 
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Fig. 31. Velocity dispersions for single stars with relative 
age errors <25% as functions of age (Fig. [201 bottom). 
From top to bottom, the total, U, V, and W velocity dis- 
persion are plotted in 10 bins with equal numbers of stars. 
The lines show fitted power laws; see text for full details. 
The youngest and oldest age bins have been excluded from 
the fits to avoid biases due to unrelaxed young structures 
and thick disk stars, respectively. 

Fig. 1301 shows no distinct jump in the velocity disper- 
sion for the oldest stars as derived by Quillen & Garnett 
I|2UU1|I from the data by Edvardsson et al. I|1993|) . Our re- 
sult, derived from a sample ~ 75 times larger and with far 
more refined procedures for the age determination, should 
clearly be more reliable. 

According to the classic Oort relation, ayjuxj should 
be constant and equal to 0.5 for a flat rotation curve (we 
find a value of - 0.63). Miihlbauer & Dehnen lj2"UU3"|l how- 
ever, show that when the true velocity dispersion and non- 
axisymmetric disturbances of the disk are taken into ac- 
count, large variations can occur. The smaller exponents 
for the in-disk heating (cy^cy) compared to the out-of- 
the-disk heating (aw) gives further constraints to be ful- 
filled by models trying to explain the observed kinematic 
heating of the disk. At least four mechanisms have been 
proposed to explain the heating: fast perturbers from the 
halo, such as 10 6 M Q black holes; slow perturbers in the 
disk, such as giant molecular clouds; large scale pertur- 
bations of the disk caused by spiral arms (De Simone et 
al. E004) or the bar (e.g. FuxEOHI}; and finally heating 
caused by infalling satellite galaxies. 

Massive black holes are improbable candidates due to 
other observational constraints, apart from the fact that 
their heating index of 0.5 is too high. Infalling satellites 
should result in a single dramatic heating of the disk, such 
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as the creation of the thick disk. Stable spiral arm patterns 
increase the random motions of stars within the plane but 
not perpendicular to it, and also become inefficient heaters 
when the epicyclic radius of a star gets larger than the 
length scale of the spiral pattern. Molecular clouds heat 
stars in all three directions, but in isolation they are in- 
efficient heaters with a total exponent of only 0.21 and a 
vertical exponent of 0.26 in the simulations of Hanninen & 
Flynn ( 2002). Clearly, further simulations including realis- 
tic descriptions of all heating components are needed, and 
indeed the recent work on the effects of stochastic, tran- 
sient spiral wave structures by De Simone et al. ( 2003) is 
able to produce exponents in the observed range. 

Widen et al. I|1996[> discussed the diffusion of stellar 
Galactic orbits as a possible explanation of the scatter 
from an otherwise single-valued age-metallicity relation, 
combined with a radial metallicity gradient in the disk. 
From our redetermination of the age-velocity relations 
based on a much larger sample and better data than avail- 
able to them, we hnd a true dispersion of velocities and 
orbits which is insufficient to support that interpretation 
(see also Nordstrom et al. I1999fl . One might suspect the 
velocity data used by Wielen et al. (19963 to be contami- 
nated by thick-disk or halo stars, or by undetected binary 
orbital motion. 

Finally, we note that the continued heating of the disk 
throughout its lifetime shown in Fig. 1311 is inconsistent 
with the results of Quillen & Garnett (2001) who found 
the heating of the thin disk to saturate after ~ 2 Gyr, 
with an abrupt increase in velocity dispersion when the 
thick disk appeared at ~ 10 Gyr. Presumably, this is due 
to their much smaller sample of stars (from Edvardsson 
et al. 11993(1 and consequently larger statistical errors than 
ours, as well as possible temperature or colour shifts of the 
evolution models used to determine ages for cool dwarfs 
(see discussion in Sect. I4.5.1|l . 

Fig. l31l does show a modest increase in velocity disper- 
sion for the oldest age bin, suggestive of the appearance 
of the thick disk. The key new result is, however, that the 
disk heating does not saturate at an early stage. Continued 
heating of the disk would appear to make it more difficult 
to trace old moving groups and other fossils of past merger 
events in the present Galactic disk (see Freeman & Bland- 
Hawthorn 12002) . 

6.5. Metallicity-velocity relations 

Given the lack of correlation between age and metallic- 
ity demonstrated in Sect. 16.21 plots of space motions vs. 
metallicity are not merely trivial transformations of the 
age- velocity relations in Sect. 16.41 but convey independent 
information. 

Fig. shows these diagrams, which display the ex- 
pected similarity with both Figs. |20] and E01 but reveal 
new, significant features. The [Fe/H] - V and [Fe/H] - U 
diagrams are especially interesting, as they show that the 
structures seen in Fig. 12*01 cover a wide range in [Fe/H] 



as well as in age (Fig. EOJ). This is quite unlike the be- 
haviour expected for stars in classical moving groups, 
which are presumed to have been born together, but re- 
sembles the kinematic structures produced by an inhomo- 
geneous Galactic potential, e.g. by stochastic spiral waves 
(De Simone et al. lSUMfl or the bar (e.g. Fux IMITj) . 

In particular, the group of stars near V = —50 km 
s^ 1 and U = —50 km s _1 and with [Fe/H] between ~ 
—0.7 and ~ +0.2 suggest the existence of a population 
of kinematically focused stars with a range of thin-disk 
metallicities in addition to the bona fide thick-disk stars. 
This possibility is of interest in studies of the detailed 
characteristics of kinematically defined samples of thick- 
disk and thin-disk stars (see Bensby et al. l2003l for a recent 
example) . 

7. Conclusions and outlook 

The results presented here are the culmination of many 
years of effort by many people. Indeed, this is in essence 
the project for which not only the Stromgren uvby photo- 
metric system fStromgren ll963(l . but also the CORAVEL 
radial-velocity scanners (Mayor 1985) were developed 
(and an important part of the science case for the 
Hipparcos satellite). We trust that the material presented 
here will remain useful to workers in the field for years to 
come. 

The properties of the data set and the results of our 
first analysis of it are described in detail above; these de- 
scriptions will not be repeated here. Instead, it is interest- 
ing to reflect on the change of scene in Galactic research 
that has taken place in the three decades since the incep- 
tion of this project (cf. the recent review by Freeman & 
Bland-Hawthorn EO^l "The New Galaxy"). 

At the time, the fundamental diagnostics of the evo- 
lution of the Galactic disk were the relations discussed 
in Sect. and deriving isochrone ages was considered 
a fairly straightforward matter. What was needed was a 
much larger and better-defined data base and an improve- 
ment in accuracy to remove what was considered to be 
merely(?) observational scatter in the relations. 

Having completed the observational effort and the first 
analysis of the data, we find ourselves with much bet- 
ter techniques for determining isochrone ages, but also 
a rather more conservative view of their accuracy than 
before. Ironically, the similar improvements in observa- 
tional accuracy (Edvardsson et al. 119931 being the stan- 
dard example) have only underscored the disagreement of 
the mean relations with the predictions of simple Galactic 
evolution models and the significance of the scatter around 
the relations in terms of the real mechanisms of the evo- 
lution of the disk. 

Our age-metallicity diagrams (Fig. I27I28|I should lay 
to rest the discussion whether closed-box models or single- 
valued predictions are adequate representations of the evo- 
lution of the Galactic disk. Models that allow for the co- 
existence of metal-poor and metal-rich stars throughout 
most of the life of the thin disk will be needed for fur- 
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ther progress. Similarly, our new age- velocity relations will 
place much stronger constraints on models of the dynam- 
ical heating of disk stars; few if any of the traditional 
candidate mechanisms seem able to match the data. The 
new results will also be of importance for assessing the 
survival rates of star clusters and merger remnants in the 
disk and perhaps allow to detect the dynamical effects of 
the bar of the Milky Way (Fux 2001). Finally, a reanalysis 
of the "G-dwarf problem" from a much larger and more 
rigorously selected data base points to a similar conclu- 
sion: The Galaxy is a far more complicated and interest- 
ing subject than ever before. The present work should lay 
the foundation for learning more about it. 

Looking ahead a decade from now, the ESA corner- 
stone mission GAIA (Perryman et al. 2001} will provide 
the next quantum leap in our knowledge of the Galaxy. 
Obtaining the complementary photometry and radial ve- 
locities needed to fully exploit the astrometric data from 
Hipparcos was left to such independent ground-based ef- 
forts as the present programme (see also Udry et al. H997f) . 
In contrast, GAIA will obtain these observations with the 
same satellite payload as the astrometric data, although 
with much larger radial- velocity errors than achieved here 
(cf. Fig. I2J) • The present material should remain useful un- 
til the results from GAIA appear, not only for studying 
our Galaxy, but also in the efforts to optimise the observ- 
ing and data reduction strategies for GAIA and for such 
precursor programmes as the RAVE survey (Steinmetz 
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Table 1. Tabic 1 will only be available in electronic 
form at the CDS. The two first pages of the table, listing 
the first 100 stars, are given at the end of this paper as 
sample of its content and format. 

Table 2. Table 2 contains the mass ratios determined for 
511 double-lined binary systems; it will be available as 
for Table 1. 
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Description of the columns in the catalogue 
(Table 1; in electronic form only) 

Left-hand pages (Fields 1-26) 

(1) HIP : Number in the Hipparcos catalogue (ESA 
I1997fl . when available. 

(2) Name : Other designation (HD, BD, CD or CP num- 
ber, in order of preference). A second number following 
a slash indicates that a double star with separate HD 
numbers has been observed together. 

(3) Comp : If the star is a member of a multiple system, 
the component (s) included in the photometry are identi- 
fied here. 

(4) f b : This flag identifies confirmed and suspected 
binaries. The information can come from one or several 
sources such as photometry, radial velocity or astrometry. 

(5) f s : This flag identifies stars in the additional sample 
of cool dwarfs south of declination —26° (see sect. 12.21 
andE3- 

(6-7) RA & Dec : Equatorial coordinates of the star 
(J2000 equinox and epoch), from the Tycho-2 catalogue. 

(8-9) I, b : Galactic coordinates of the star. 

(10) V : Johnson V magnitude from the published uvby 
photometry. 

(11) b-y : The b-y colour in the Stromgren uvby system. 

(12) 13 : The Hp index in the Stromgren system. 

(13) E(b-y) : Colour excess, from the calibrations of 
Olsen ltT9m 

(14) \ogT e ff : Logarithm of the effective temperature, 
from the uvby indices and the calibration by Alonso et 
al. (JUSl). 

(15) [Fe/H] : Metallicity of the star as determined from 
the uvby(3 indices and the calibrations by Schuster & 
Nissen (:1982J), Edvardsson et al. 1)1998(1 or the new 
calibration defined in Sect. 14.31 

(16) d : Distance of the star, in parsecs (pc). Sect. 14.41 
explains why some stars have no entry in this column. 

(17) My ■ Absolute magnitude, from the apparent 
magnitude, distance, and colour excess. 

(18) SMy '. The magnitude difference between the star 
and the ZAMS. 

(19) f r : Source for the distance (see text). H = Hipparcos 
parallax; F,G = F or G-star photometric distance. 

(20) f g : Suspected giant flag. An asterisk (*) indicates 
a disagreement between the photometric determination 
and the Hipparcos parallax at the 3cr level, suggesting 
that the star is a giant not detected from the photometry. 

(21-23) age, a l ™ e , af g f : The age determined for the 
star, if any, with upper and lower la confidence limits, 
all in Gigayears. 

(24-26) mass, a l ^ ss , a^ s : The mass determined for 
the star, with upper and lower ler confidence limits, all in 
Solar masses. 



Right-hand pages (Fields 27-51) 

(27) Name : Repeated from field 2. 

(28) V r : Mean radial velocity of the star. For double- 
lined binaries, the computed systemic velocity is given if 
so indicated by the flag fd in field 34. 

(29) m.e. : Mean error of the mean radial velocity. See 
text for details. 

(30) <7y r : Standard deviation of the individual radial- 
velocity measurements if TV > 1. 

(31) TV : The number of individual radial- velocity 
measurements. 

(32) AT : Time-span, in days, from the first to the last 
radial-velocity observation. 

(33) P(x 2 ) ■ Probability that the observed scatter of the 
radial velocities is due to random observational errors 
only. 

(34) fd : This flag indicates that the star is a spectro- 
scopic binary with a systemic radial velocity given in field 
28. 

(35) f v : Source of the radial velocity: C = Coravel; A = 
CfA; L = literature. 

(36) vsmi : Rotational velocity of the star, in km s _1 . 

(37-39) jj, a * /is c M : Proper motion in right ascension and 
declination and standard error of the combined motion, 
in milliarcsec/year. Mostly from the Tycho-2 catalogue. 

(40-41) 7r ov : Hipparcos parallax and its standard error, 
in milliarcsec. 

(42-44) U V W : Heliocentric space velocity components 
for the star (U positive towards the Galactic centre), in 
km s _1 . 

(45-46) Rgai , z : Galactic position of the star, in kpc. 
(47-50) R m in Rmax e z max : Perigalactic and apogalactic 
distance, eccentricity, and maximum distance from the 
galactic plane of the computed galactic orbit, in kpc. 

(51) /„ : General note (see below). 

a: Double star with Am < 5 mag 

b: Double star with Am < 5 mag from Hipparcos 

c: Variable star 

d: Simbad note 

e: See Olsen ifTflT^ 

f: See Olsen lETTO 

g: See Olsen (|T3gO) 

h: See Olsen itTMfl 

i: See Olsen (|1994aJ) 

j: See Olsen l|1994hll 

k: See Abt et al. lfT9T9")l 

1: See Abt ipgl) 

m: See Abt ifT^SB) 

n: See Gray & Garrison l)1989|l 

o: See Gray fTT^ 

p: See Henry et al. ifTMfll 
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Table 2. Sample left-hand page of the catalogue (Fields 1-25 for the first 100 stars). 
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Table 2. Sample right-hand page of the catalogue (Fields 26-50 for the first 100 stars). 
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